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1 Introduction 

One of the main open problems in astrophysics is the dominance of dark en- 
tities, the dark matter and the dark energy, in the existing description of the 
universe given by the standard ACDM cosmological model [1] based on the 
cosmological principle (homogeneity and isotropy of the space-time), which 
selects the class of Friedmann-Robertson- Walker (FWR) space-times. After 
the transition from quantum cosmology to classical astrophysics, with the 
Hcisenberg cut roughly located at a suitable cosmic time (ss 10 5 years after 
the big bang) and at the recombination surface identified by the cosmic mi- 
crowave background (CMB), one has a description of the universe in which 
the known forms of baryonic matter and radiation contribute only with a few 
percents of the global budget. One has a great variety of models trying to ex- 
plain the composition of the universe in accelerated expansion (based on data 
on high red-shift supernovae, galaxy clusters and CMB): WIMPS (mainly 
super-symmetric particles), f(R) modifications of Einstein gravity (with a 
modified Newton potential), MOND (with a modification of Newton law),... 
for dark matter; cosmological constant, string theory, back-reaction (spatial 
averages, non-linearity of Einstein equations), inhomogeneous space-times 
(Lemaitre-Tolman-Bondi, Szckeres), scalar fields (quintessence, k-essence, 
phantom), fluids (Chaplygin fluid), .... for dark energy. 

Most of these developments rest on a description based on a family of 
FRW space-times with nearly flat 3-spaces (as required by CMB data) as the 
reference space-times where to interpret the astronomical data (luminosity, 
light spectrum, angles) on the 2-dimensional sky vault. Therefore, the starting 
point is the extension of the standards of relativistic metrology near the Earth 
and in the Solar System to astronomy: to reconstruct a 4-dimensional space- 
time one needs new standards of time and length like the cosmic time and the 
luminosity distance (or any other astrometric definition, see Ref. [2] ) allowing 
to define an International Celestial Reference System (ICRS) [3], namely 
a 4-coordinate system describing a 3-universe evolving in time, where the 
astronomical data have to be dynamically interpreted according to Einstein 
gravity or some of its extensions. 

The aim of this Lecture is to suggest a new viewpoint on the origin of 
dark matter, and maybe also of dark energy, starting from a re-reading of 
the general covariance of Einstein general relativity (GR), which could be 
also applied to every generally covariant extension of this theory if needed. It 
is an extended version of the review paper [4]. In this Introduction one will 
delineate the framework of our approach and then in the subsequent Sections 
one will give more details of the various topics. 

The gauge group of the Lagrangian formulation of Einstein GR, the diffeo- 
morphism group, implies that the 4-coordinates of the space-time are gauge 
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variables. As a consequence, the search of GR observables is restricted to 4~ 
scalars and at the theoretical level one tries to describe gravitational dynam- 
ical properties in term of them. However, inside the Solar System the exper- 
imental localization of macroscopic classical objects is unavoidably done by 
choosing some convention for the local 4-coordinates of space-time. Atomic 
physicists, NASA engineers and astronomers have chosen a series of refer- 
ence frames and standards of time and length suitable for the existing tech- 
nology [5]. These conventions determine certain Post-Minkowskian (PM) 4- 
coordinate systems of an asymptotically Minkowskian space-time, in which 
the instantaneous 3-spaces are not strictly Euclidean. Then these reference 
frames are seen as a local approximation of a reference frame in ICRS, where 
however the space-time has become a cosmological FWR one, which is only 
conformally asymptotically Minkowskian at spatial infinity. A search of a 
consistent patching of the 4-coordinates from inside the Solar System to the 
rest of the universe will start when the data from the future GAIA mission 
[6] for the cartography of the Milky Way will be available. This will allow a 
PM definition of a Galactic Reference System containing at leat our galaxy. 
Let us remark that notwithstanding the FRW instantaneous 3-spaces are not 
strictly Euclidean, all the books on galaxy dynamics describe the galaxies by 
means of Kepler theory in Galilei space-time. 

This state of affairs requires to revisit Einstein GR to see whether it is 
possible to identify which components of the 4-metric tensor are connected 
with the gauge freedom in the choice of the 4-coordinates and which ones de- 
scribe the dynamical degrees of freedom of the gravitational field. Since this 
cannot be done at the Lagrangian level, one must restrict himself to the class 
of globally hyperbolic, asymptotically flat space-times allowing a Hamilto- 
nian description starting from the description of Einstein GR in terms of the 
ADM action [7] instead than in terms of the Einstein-Hilbert one. In canon- 
ical ADM gravity one can use Dirac theory of constraints [8] to describe 
the Hamiltonian gauge group, whose generators are the first-class constraints 
of the model. The basic tool of this approach is the possibility to find so- 
called Shanmugadhasan canonical transformations [9] , which identify special 
canonical bases adapted to the first-class constraints (and also to the second- 
class ones when present). In these special canonical bases the vanishing of 
certain momenta (or of certain configurational coordinates) corresponds to 
the vanishing of well defined Abelianized combinations of the first-class con- 
straints (Abelianized because the new constraints have exactly zero Poisson 
brackets even if the original constraints were not in strong involution). As 
a consequence, the variables conjugate to these Abelianized constraints are 
inertial Hamiltonian gauge variables describing the Hamiltonian gauge free- 
dom. The remaining 2+2 conjugate variables describe the dynamical tidal 
degrees of freedom of the gravitational field (the two polarizations of gravi- 
tational waves in the linearized theory). If one would be able to include all 
the constraints in the Shanmugadhasan canonical basis, these 2+2 variables 
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would be the Dirac observables of the gravitational field, invariant under the 
Hamiltonian gauge transformations. However such Dirac observables are not 
known: one only has statements bout their existence [10]. Moreover, in gen- 
eral they are not 4-scalar observables. The problem of the connection between 
the 4-diffeomorphism group and the Hamiltonian gauge group was studied in 
Ref.[ll] by means of the inverse Legendre transformation and of the notion 
of dynamical symmetry. The conclusion is that on the space of solutions of 
Einstein equations there is an overlap of the two types of observables: there 
should exists special Shanmugadhasan canonical bases in which the 2+2 Dirac 
observables become 4-scalars when restricted to the space of solutions of the 
Einstein equations. In any case the identification of the inertial gauge com- 
ponents of the 4-metric is what is needed to make a fixation of 4-coordinates 
as required by relativistic metrology. 

Another problem is that asymptotically flat space-times have the SPI 
group of asymptotic symmetries (direction-dependent asymptotic Killing 
symmetries) [12] and this is an obstruction to the existence of asymptotic 
Lorentz generators for the gravitational field [13]. However if one restricts the 
class of space-times to those not containing super-translations [14], then the 
SPI group reduces to the asymptotic ADM Poincare group [15]: these space- 
times are asymptotically Minkowskian, they contain an asymptotic Minkowski 
4-metric (to be used as an asymptotic background at spatial infinity in the 
linearization of the theory) and they have asymptotic inertial observers at 
spatial infinity whose spatial axes may be identified by means of the fixed 
stars of star catalogues 1 . Moreover, in the limit of vanishing Newton constant 
(G = 0) the asymptotic ADM Poincare generators become the generators of 
the special relativistic Poincare group describing the matter present in the 
space-time. This is an important condition for the inclusion into GR of the 
classical version of the standard model of particle physics, whose properties 
are all connected with the representations of this group in the inertial frames 
of Minkowski space-time. In absence of matter a sub-class of these space- 
times is the (singularity-free) family of Chrstodoulou-Klaincrmann solutions 
of Einstein equations [16] (they are near to Minkowski space-time in a norm 
sense and contain gravitational waves). 

Moreover, in this restricted class of space-times the canonical Hamiltonian 
is the ADM energy [17], so that there is no frozen picture like in the "spatially 
compact space-times without boundaries" used in loop quantum gravity 2 . 



1 The fixed stars can be considered as an empirical definition of spatial infinity of the 
observable universe. 

2 In these space-times the canonical Hamiltonian vanishes and the Dirac Hamiltonian is a 
combination of first-class constraints, so that it only generates Hamiltonian gauge transfor- 
mations. In the reduced phase space, quotient with respect the Hamilonian gauge group, 
the reduced Hamiltonian is zero and one has a frozen picture of dynamics. This class of 
space-times fits well with Machian ideas (no boundary conditions) and with interpretations 
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To take into account the fermion fields present in the standard particle 
model one must extend ADM gravity to ADM tetrad gravity . Since our 
class of space-times admits orthonormal tetrads and a spinor structure [19], 
the extension can be done by simply replacing the 4-metric in the ADM 
action with its expression in terms of tetrad fields, considered as the basic 16 
configurational variables substituting the 10 metric fields. 

To study ADM tetrad gravity the preliminary problem is to choose a co- 
ordinatization of the space-time compatible with relativistic metrology. This 
requires a definition of global non-inertial frames, because the equivalence 
principle forbids the existence of global inertial frames in GR. Due to the 
Lorentz signature of the space-time this is a non-trivial task already in spe- 
cial relativity (SR): there is no notion of instantaneous 3-space, because the 
only intrinsic structure is the conformal one, i.e. the light-cone as the locus of 
incoming and outgoing radiation. The existing coordinatizations, like either 
Fermi or Riemann-normal coordinates, hold only locally They are based on 
the 1+3 point of view, in which only the world-line of a time- like observer is 
given. In each point of the world-line the observer 4-velocity determines an 
orthogonal 3-dimcnsional space-like tangent hyper-plane, which is identified 
with an instantaneous 3-space. However, these tangent planes intersect at 
a certain distance from the world-line (the so-called acceleration length de- 
pending upon the 4- acceleration of the observer [20]), where 4-coordinates of 
the Fermi type develop a coordinate singularity. Another type of coordinate 
singularity is developed in rigidly rotating coordinate systems at a distance 
r from the rotation axis where co r = c (uj is the angular velocity and c the 
two-way velocity of light). This is the so-called "horizon problem of the ro- 
tating disk" : a time- like 4-velocity becomes a null vector at wr = c, like it 
happens on the horizon of a black-hole. See Ref.[21] for a classification of the 
possible pathologies of non-inertial frames and on how to avoid them. 

In this Lecture one will review the way out from these problems based on 
the 3+1 point of view in which, besides the world-line of a time-like observer, 
one gives a global nice foliation of the space-time with instantaneous 3-spaces. 
Then a metrology-oriented notion of 4-coordinates, the so-called radar 4~ 
coordinates first introduced by Bondi [22], is introduced in these global non- 
inertial frames. One will give the conditions for a foliation to be nice, i.e. for 
the absence of pathologies like the ones of the rotating disk and of the Fermi 
coordinates. 

Let us remark that the theory of global non-inertial frames is also needed 
to speak of predictability in a (either classical or quantum) theory in which 
the basic equations of motion are partial differential equations (PDE). To be 
able to use the existence and unicity theorem for the solutions of PDE's, one 



in which there is no physical time like the one in Ref. [18]. However, it is not clear how to 
include in this framework the standard model of particle physics. 
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needs a well-posed Cauchy problem, whose prerequisite is a sound definition 
of an instantaneous 3-space (i.e. of a clock synchronization convention) where 
the Cauchy data are given. To give the data on a space-like surface is not 
factual, but with the data on the backward light-cone of an observer it is 
not yet possible to demonstrate the theorem. However, also the 1+3 point of 
view is non factual, because it requires the knowledge of a world-line from 
the whole past to all the future. 

A Section of this Lecture will be devoted to the developments in relativistic 
particle mechanics made possible by the 3+1 point of view in SR [23], [21], 
[24]. By means of parametrized Minkowski theories [23], [21], one can get 
the description of arbitrary isolated systems (particles, strings, fluids, holds) 
admitting a Lagrangian formulation in arbitrary non-inertial frames with the 
transition among non-inertial frames described as a "gauge transformation" 
(general covariance under the frame-preserving diffeomorphisms of Ref. [25]). 
Moreover this framework allows us to define the inertia! and non-inertial 
rest frames of the isolated systems, where to develop the rest-frame instant 
form of the dynamics and to build the explicit form of the Lorentz boosts 
for interacting systems. This makes possible to study the problem of the 
relativistic center of mass [26], relativistic bound states [27, 28, 29], relativistic 
kinetic theory and relativistic micro-canonical ensemble [30] and various other 
systems [31, 32]. Moreover a Wigner-covariant relativistic quantum mechanics 
[33], with a solution of all the known problems introduced by SR, has been 
developed after some preliminary work done in Ref. [34] . This will allow us to 
study relativistic entanglement. 

After this digression in SR one defines global non-inertial frames with 
radar 4-coordinates in the asymptotically Minkowskian space-times of GR 3 
and one gives the parametrization of the tetrads and of the 4-mctric in them. 
The absence of super-translations implies that these non-inertial frames are 
non-inertial rest frames of the 3-universc. Starting from the ADM action 
for tetrad gravity one defines the Hamiltonian formalism in a phase space 
containing 16 configurational field variables and 16 conjugate moments. One 
identifies the 14 first-class constraints of the system and one finds that the 
canonical Hamiltonian is the weak ADM energy (it is given as a volume in- 
tegral over 3-space). The existence of these 14 first-class constraints implies 
that 14 components of the tetrads (or of the conjugate momenta) are Hamil- 
tonian gauge variables describing the inertia! aspects of the gravitational field 
(6 of these inertial variables describe the extra gauge freedom in the choice 
of the tetrads and in their transport along world- lines) . Therefore there are 
only 2+2 degrees of freedom for the description of the tidal dynamical as- 
pects of the gravitational field. The asymptotic ADM Poincare generators 

3 While in SR Minkowski space-time is an absolute notion, unifying the absolute notions 
of time and 3-space of the non-rclativistic Galilei space-time, in GR there is no absolute 
notion: space-time becomes dynamical [35] with its metric structure satisfying Einstein 
equations. 



8 



Luca Lusanna 



can be evaluated explicitly. Till now the type of matter studied in this frame- 
work consists of the electro-magnetic field and of N charged scalar particles, 
whose signs of the energy and electric charges are Grassmann- valued to reg- 
ularize both the gravitational and electro-magnetic self-energies (it is both a 
ultraviolet and an infrared regularization), 

Then it will be shown that there is a Shanmugadhasan canonical transfor- 
mation [36] (implementing the so-called York map [37] and diagonalizing the 
York-Lichnerowics approach [38]) to a so-called York canonical basis adapted 
to 10 of the 14 first-class constraints. Only the super-Hamiltonian and super- 
momentum constraints, whose general solution is not known, are not included 
in the basis, but it is clarified which variables are to be determined by their 
solution. Among the inertial gauge variables there is the York time [39] 3 K , 
i.e. the trace of the extrinsic curvature of the 3-spaces as 3-manifolds embed- 
ded into the space-time. It is the only gauge variable which is a momentum in 
the York canonical basis 4 : this is due to the Lorentz signature of space-time, 
because the York time and three other inertial gauge variables can be used 
as 4-coordinates of the space-time (see Ref. [35] for this topic and for its rele- 
vance in the solution of the hole argument). Therefore an identification of the 
inertial gauge variables to be fixed to get a 4-coordinate system in relativistic 
metrology was found. In the first paper of Ref. [40] there is the expression of 
the Hamilton equations for all the variables of the York canonical basis. 

An important remark is that in the framework of the York canonical basis 
the natural family of gauges is not the harmonic one, but the family of 3- 
orthogonal Schwinger time gauges in which the 3-metric in the 3-spaces is 
diagonal. 

Both in SR and GR an admissible 3+1 splitting of space-time has two 
associated congruences of time-like observers [21], geometrically defined and 
not to be confused with the congruence of the world-lines of fluid elements, 
when relativistic fluids are added as matter in GR [41, 42, 43]. One of the 
two congruences, with zero vorticity, is the congruence of the Eulerian ob- 
servers, whose 4-velocity field is the field of unit normals to the 3-spaces. 
This congruence allows us to re-express the non- vanishing momenta of the 
York canonical basis in terms of the expansion (9 = — K) and of the shear 
of the Eulerian observers. This allows us to compare the Hamilton equa- 
tions of ADM canonical gravity with the usual first-order non-Hamiltonian 
ADM equations deducible from Einstein equations given a 3+1 splitting of 
space-time but without using the Hamiltonian formalism. As a consequence, 
one can extend our Hamiltonian identification of the inertial and tidal vari- 
ables of the gravitational field to the Lagrangian framework and use it in 
the cosmological (conformally asymptotically flat) space-times: in them it is 
not possible to formulate the Hamiltonian formalism but the standard ADM 



4 Instead in Yang-Mills theory all the gauge variables are configurational. 
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equations are well denned. The time inertial gauge variable needed for rela- 
tivistic metrology is now the expansion of the Eulcrian observers of the given 
3+1 splitting of the globally hyperbolic cosmological space-time. 

The next step (see the second paper of Rcf. [40]) is the definition of a PM 
linearization of ADM tetrad gravity in the family of 3-orthogonal Schwinger 
time gauges in which one chooses 3-coordinates diagonalizing the 3-metric in 
the 3-spaces and an arbitrarily given numerical function for the York time 3 K . 
The cosmic time T cosm has to be chosen so that the 3-spaces T cosm = const. 
have an extrinsic curvature with the given value of 3 K. This PM linearization 
uses the asymptotic Minkowski 4-metric as an asymptotic background, so that 
one never splits the 4-mctric with respect to a fixed Minkowski metric in the 
bulk like in the standard approach to gravitational waves. A ultraviolet cutoff 
on the matter is needed. 

This leads to a PM formulation of gravitational waves in non-harmonic 
3-orthogonal gauges. All the constraints can be solved, an explicit expres- 
sion of the PM 4-metric can be given and the explicit form of the Hamilton 
equations for the tidal degrees of freedom of the gravitational field and the 
matter can be obtained. It is non-trivial to show that all the standard re- 
sults about gravitational waves in harmonic gauges [44] can be reproduced 
in the 3-orthogonal gauges with the help of the formalism of Ref.[45]. As 
shown in the third paper of Rcf. [40] (where the matter is restricted only to 
scalar particles), all the 4- and 3- curvature tensors of the space-time can 
be explicitly evaluated and the time-like and null geodesies can be studied. 
It is also possible to evaluate the red-shift of light rays and the luminosity 
distance finding their dependence on the York time and verify the old Hubble 
red-shift-distance law (see ref.[46]), which becomes the usual Hubble law (a 
velocity-distance relation) when one uses the standard cosmological model. 
In the Solar System the results in the 3-orthogonal gauges are compatible 
with the ones in the harmonic gauges used in relativistic metrology [5] . 

The main important result or this lecture are the PM Hamilton equations 
and the implied PM second-order equations of motion for the particles. Their 
PN limit identifies the Newton forces acting on the particles at the lowest or- 
der augmented with 1PN forces compatible with the known results on binaries 
in harmonic gauges [47]. However, there are extra 0.5PN forces, depending 
linearly on the non-local York time 3 /C = 3 K (A is the asymptotic Lapla- 
cian of the 3-space), representing either a friction or an anti-friction force 
according to the sign of 3 /C. These 0.5PN forces are our main result, because 
their effect can be re-intcrprctcd as an extra effective (time-, position- and 
velocity-dependent) contribution to the inertial mass of the particles in the 
equations used in the three main signatures for the existence of dark matter: 
the rotation curves of spiral galaxies [48] and the masses of clusters of galaxies 
from the virial theorem [49] and from weak gravitational lensing [50] , [49] . 
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While gravitational and inertial masses are equal in Einstein GR, the PM 
limit, followed by the PN one, shows that the non-Euclidean nature of the 
3-spaces implies a breaking of the Newtonian equality of the two types of 
masses, which holds only in the absolute Euclidean 3-space of Galilei space- 
time. 

As a consequence the data on dark matter can be re-read as a partial 
fixation of the non-local York time 3 IC. However to fix the York time 3 K = 
A 3 /C one needs a global information on 3 /C on the whole 3-space, in particular 
in the voids among galaxy clusters. 

Therefore one has an indication that (at least part of) dark matter could 
be re- absorbed in a PM extension of the conventions in the existing ICRS, 
such that the 3-spaces tjcrs — const, determined by a suitable ICRS time 
have a York time 3 K such that the derived non-local York time reproduces 
the data for the signatures of dark matter. 

In the Conclusions it will be suggested that also the open problem of dark 
energy could be rephrased as the determination of a suitable York time in 
inhomogeneous cosmological space-times. Therefore there is the possibility of 
an understanding of the " dark" aspects of the universe in terms of relativistic 
metrology. 
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2 Relativistic Metrology 

As shown in Ref. [51] modern relativistic metrology is not only deeply rooted 
in Maxwell theory and its quantization but is also beginning to take into 
account GR. 

The basic metrological conventions on the Earth surface are: 

a) An atomic clock as a standard of time. The fundamental conceptual time 
scale is the SI atomic second: it is the duration of 9 192 631 770 periods of the 
radiation corresponding to the transition between the two hyperhne levels of 
the ground state of the cesium 133 atom. This definition refers to a cesium 
atom at rest at a temperature of 0°. However the practical time standard is 
the International Atomic Time (TAI), which is defined as a suitable weighted 
average of the SI time kept by over 200 atomic clocks in about 70 national 
laboratories worldwide. Time scales connected with TAI are the GPS Station 
Time and the Universal Time (UC) based on Earth's rotation 5 . All the other 
existing time scales inside the Solar System are connected to this standard 
by fixed conventions. 

b) The 2- way velocity of light (only one clock is involved in its definition) , 
fixed to the value c = 299 792 458m s -1 , in place of the standard of length 
6 . To measure the 3-distance between two objects in an inertial frame one 
sends a ray of light from the first object, to which is associated an atomic 
clock, to the second one, where it is reflected and then reabsorbed by the 
first object. The measure of the flight time and the 2-way velocity of light 
determine the 3-distance between the objects. This convention is compatible 
with the Euclidean 3-space of inertial frames in Minkowski space-time. When 
the technology will allow one to measure the deviations from Euclidean 3- 
space implied by PN gravity one will need a modified convention taking into 
account a general relativistic notion of length. 

Given these standards one can think to the Global Positioning System 
(GPS) as a local standard of space-time. To define GPS one needs a conven- 
tional reference frame centered on a given time-like observer. Inside the Solar 
System one has well defined conventions for the following reference frames: 

A) The description of satellites around the Earth is done by means of 
NASA coordinates either in the International Terrestrial Reference System 
(ITRS; it is a frame fixed on the Earth surface) or in the Geocentric Celestial 
Reference System (GCRS) centered on the world-line of the Earth center (see 



5 It is based on Very Long Baseline Interferometry (VLBI) observations of distant quasars, 
on Lunar Laser Ranging (LLR) and on determination of GPS satellite orbits. 

6 The meter is the length of the path traveled by light in vacuum during a time interval 
of 1/c of a second. 
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Ref.[5]). Both of them use a geocentric coordinate time t G connected to TAI. 

B) The description of planets and other objects in the Solar System uses 
the Barycentric Celestial Reference System (BCRS), centered in the barycen- 
ter of the Solar System (see Ref.[5]). It uses a barycentric coordinate time ts 
connected to to and TAI. 

While ITRS is essentially realized as a non-relativistic non-inertial frame 
in Galilei space-time, BCRS is defined as a quasi- inertial frame, non-rotating 
with respect to some selected fixed stars, in Minkowski space-time with 
nearly-Euclidean 3-spaces (one ignores the perturbations induced from the 
Milky Way). It can also be considered as a PM space-time with 3-spaces 
having a very small extrinsic curvature of order c -2 . GCRS is obtained from 
BCRS by taking into account Earth's rotation around the Sun with a suitable 
Lorentz boost with corrections from PN gravity 7 . By taking into account the 
extension of the geoid and Earth revolution around its axis one goes from 
the quasi-Minkowskian GCRS to the quasi-Galilean ITRS. 

New problems emerge by going outside the Solar System. In astronomy 
the positions of stars and galaxies are determined from the data (luminosity, 
light spectrum, angles) on the sky, i.e. on a 2-dimcnsional spherical surface 
around the Earth with the relations between it and the observatory on the 
Earth done with GPS. To get a description of stars and galaxies as living 
in a 4-dimensional space-time one introduces the International Celestial Ref- 
erence System ICRS (see Refs. [3]). Its time scale is a "second" connected 
to GPS, TAI and SI and therefore to t G and ts- ICRS has the origin in 
the solar system barycenter, which is considered as quasi-inertial observer 
carrying a quasi-inertial (essentially non-relativistic) reference frame with 
rectangular 3-coordinates in a nearly Galilei space-time whose 3-spaces are 
nearly Euclidean. The directions of the spatial axes are effectively defined 
by the adopted coordinates of 212 extragalactic radio sources observed by 
VLBI . These radio sources (quasars and AGN, active galactic nuclei) are 
assumed to have no observable intrinsic angular momentum. Thus, the ICRS 
is a space-fixed system, more precisely a kinematically non-rotating system, 
which provides the orientation of BCRS. 

In astronomy the unit of length is the astronomical unit A U, approximately 
equal to the mean Earth-Sun distance. Measurements of the relative positions 
of planets in the Solar System are done by radar: one measures the time taken 
for light to be reflected from an object using the conventional value of the 
velocity of light c. Both for objects inside the Solar System and for the nearest 
stars one measure the distance with the trigonometric parallax by using the 
propagation of light and its velocity c in inertial frames. One measures the 
difference (the inclination angle) in the apparent position of an object viewed 



7 See Ref. [52] for possible gravitational anomalies inside the Solar System. 
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along two different lines of sight at two different times and then uses Euclidean 
geometry to evaluate the distance. The used unit in astronomy is the parsec, 
which is 3.26 light-years or 3.26 10 16 meters. 

This convention cannot be used for more distant either galactic or extra- 
galactic objects. New notions like standard candles, dynamical parallax, spec- 
troscopic parallax, luminosity distance, are needed [2]. These notions in- 
volve both aspects of light propagation in curved space-times and cosmolog- 
ical assumptions like the Hubble law (velocity- red- shift linear relation). 

However if one takes into account the description of the universe given 
by cosmology, the actual cosmological space-time cannot be a nearly Galilei 
space-time but it must be a cosmological space-time with some theoretical 
cosmic time. In the standard cosmological model [1] it is a homogeneous and 
isotropic FRW space-time whose instantaneous 3-spaces have nearly van- 
ishing internal 3-curvature, so that they may locally be replaced with Eu- 
clidean 3-spaces as it is done in galactic dynamics. However they have a 
time-dependent conformal factor (it is one in Galilei space-time) responsi- 
ble for the Hubble constant regulating the expansion of the universe. As a 
consequence the transition from the astronomical ICRS to an astrophysical 
description taking into account cosmology is far from being understood. 

What is still lacking is a PM extension of the celestial frame such that the 
PM BCRS frame is its restriction to the solar system inside our galaxy. In 
particular one needs the definition of a coordinate time ticns connected to ts 
such that the 3-spaces ticns = const, have a very small internal 3-curvature 
and a suitable extrinsic curvature as sub-manifolds of the space-time con- 
nected with the Hubble constant. In this way this astronomical PM ICRS 
would be consistent with the FRW cosmological space-times used in astro- 
physics except for the conformal factor determining the accelerated expansion 
of the universe and creating problems in the metrological use of fixed stars. 

Hopefully at least an PM extension of ICRS including our galaxy (with 
the definition of a galactic coordinate system) will be achieved with the ESA 
GAIA mission devoted to the cartography of the Milky Way [6] . 
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3 Clock Synchronization and Global Non-Inertial 
Frames in Minkowski Space-Time 

Since in the Minkowski space-time of SR time is not absolute, there is no 
intrinsic notion of 3-space and of synchronization of clocks: both of them have 
to be defined with some convention. As a consequence the 1-way velocity of 
light from one observer A to an observer B has a meaning only after a choice 
of a convention for synchronizing the clock in A with the one in B. Therefore 
the crucial quantity in special relativity is the 2-way (or round trip) velocity 
of light c involving only one clock. It is this velocity (a kind of mean velocity) 
which is isotropic and constant in SR and replaces the standard of length in 
relativistic metrology. 

Einstein convention for the synchronization of clocks in Minkowski space- 
time uses the 2-way velocity of light to identify the Euclidean 3-spaces of the 
inertial frames centered on an inertial observer A by means of only its clock. 
The inertial observer A sends a ray of light at x° towards the (in general 
accelerated) observer B; the ray is reflected towards A at a point P of B 
world-line and then reabsorbed by A at x°j-, by convention P is synchronous 
with the mid-point between emission and absorption on A's world-line, i.e. 
x° P = x° + 2 (Xf — x°) = \ (x° + x°f). This convention selects the Euclidean 
instantaneous 3-spaces x° = ct = const, of the inertial frames centered on 
A. Only in this case the one-way velocity of light between A and B coincides 
with the two-way one, c. However, as said in the Introduction, if the observer 
A is accelerated, the convention can breaks down due the possible appearance 
of coordinate singularities. 

As a consequence, a theory of global non-inertial frames in Minkowski 
space-time has to be developed in a metrology-oriented way to overcame the 
pathologies of the 1+3 point of view. This has been done in the papers of 
Ref.[21] based on the 3+1 point of view and on the use of observer-dependent 
Lorentz scalar radar 4-coordinates. This theory and its implications for the 
description of isolated systems in SR will be reviewed in this Section. 



3.1 3+1 Splittings of Minkowski Spacetime and Radar 
4- Coordinates 

Assume that the world- line a+ (r) of an arbitrary time-like observer carrying 
a standard atomic clock is given: t is an arbitrary monotonically increasing 
function of the proper time of this clock. Then one gives an admissible 3+1 
splitting of Minkowski space-time, namely a nice foliation with space-like 
instantaneous 3-spaces S T . It is the mathematical idealization of a protocol 
for clock synchronization: all the clocks in the points of E T sign the same 
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time of the atomic clock of the observer 8 . On each 3-space S T one chooses 
curvilinear 3-coordinates o~ r having the observer as origin. These are the 
Lorentz-scalar and observer-dependent radar 4~ coordinates a A = (t;<j t ). 

If x^ h- > a A (x) is the coordinate transformation from the Cartesian 4- 
coordinates x^ of a reference inertial observer to radar coordinates, its inverse 
a A i y x^ = z^(t, a r ) defines the embedding functions z^(t, a r ) describing the 

3- spaces S T as embedded 3-manifold into Minkowski space-time. The induced 

4- metric on S T is the following functional of the embedding 4 gAB{i~,o- r ) = 

[z^r]^ z^](t, a r ), where z\ = dz^/da A and 4 ?7 Aily = e (H ) is the 

flat metric 9 . While the 4- vectors z^(r,a u ) are tangent to E T , so that 
the unit normal l^{r,a u ) is proportional to e A ' Q( 9 7 [zf z% zJ](t, a u ), one has 
z£(r,CT r ) = [Nl»+N r z?}(T,o- r ) withiV(r,a r ) = e [z? Z m ](t, a r ) = l+n( T ,a r ) 
and N r (r 1 o- r ) = — eg rr (r, er r ) being the lapse and shift functions. 

As a consequence, the components of the 4-metric ^^(t, er r ) have the 
following expression 

e 4 5T r = N 2 - N r N r , -e 4 5rr = N r = 3 g rs N s , 

3 

3 4 x 3 3 

Qrs £ 9rs / J e (a)r e (a)s 

a=l 

3 

= 4> 2/3 E e ' Ei=1 7Sa Ri v ™( 01 ) ^ l )> (!) 

a=l 

where 3 e( a ) r (r, a") are cotriads on ZY, 4' 2 {T 1 o- r ) — det 3 g rs (T,o- r ) is the 3- 
volume element on Z r , A a (r, a r ) — [(j) 1 / 3 e^s=i 76 " Ri ](t, a r ) are the posi- 
tive eigenvalues of the 3-metric (j aa are suitable numerical constants) and 
V(9 1 (t, cr r )) are diagonalizing rotation matrices depending on three Euler 
angles. 

Therefore starting from the four independent embedding functions z^ (r, a r ) 
one obtains the ten components 4 gAB of the 4-metric (or the quantities N, 
N r , <f>, R a , 9 l ), which play the role of the inertial potentials generating the 
relativistic apparent forces in the non- inertial frame. It can be shown [21] 
that the usual non-relativistic Newtonian inertial potentials are hidden in 
the functions N, N r and 9 l . The extrinsic curvature tensor 3 K rs (r, a u ) = 
[ {N r \ s + N s \ r — d T 3 g rs )](T, a u ), describing the shape of the instantaneous 
3-spaccs of the non-inertial frame as embedded 3-sub-manifolds of Minkowski 



8 It is the non-factual idealization required by the Cauchy problem generalizing the existing 
protocols for building coordinate system inside the future light-cone of a time-like observer. 

9 e = ±1 according to either the particle physics e = 1 or the general relativity e = — 1 
convention. 
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space-time, is a secondary inertial potential, functional of the ten inertial po- 
tentials 4 gAB- 

The foliation is nice and admissible if it satisfies the conditions: 

1) N(t, <r r ) > in every point of S T so that the 3-spaces never intersect, 
avoiding the coordinate singularity of Fermi coordinates; 

2) e 4 g TT (T, <r r ) > 0, so to avoid the coordinate singularity of the rotat- 
ing disk, and with the positive-definite 3-metric 3 g rs (r, a u ) = —e i g rs (T,a u ) 
having three positive eigenvalues (these are the M0ller conditions [53]); 

3) all the 3-spaces S T must tend to the same space-like hyper-plane at 
spatial infinity with a unit normal e^, which is the time- like 4- vector of a set of 
asymptotic ortho-normal tetrads e^. These tetrads are carried by asymptotic 
inertial observers and the spatial axes e£f are identified by the fixed stars of 
star catalogues. At spatial infinity the lapse function tends to 1 and the shift 
functions vanish. 



3.2 Global Non-Inertial Frames in Minkowski 
Spacetime 

By using the asymptotic tetrads one can give the following parametrization 
of the embedding functions 

z"{t, a r ) = x»{t) + e\ F a (t, a r ), F a (t, 0) = 0, 

x»(T)=xt + e>Xf A (T), (2) 

where x^(t) is the world- line of the observer. The functions / (r) deter- 
mine the 4- velocity u^(t) = x^ (t) / e x 2 (t) (x m (t) = ) and the 4- 
acceleration a Ai (r) = d% 1 ^ of the observer. 

The M0ller conditions are non-linear differential conditions on the func- 
tions f A (r) and F a (t, <r r ), so that it is very difficult to construct explicit 
examples of admissible 3+1 splittings. When these conditions are satisfied 
Eqs.(2) describe a global non- inertial frame in Minkowski space-time. 

Till now the solution of M0ller conditions is known in the following two 
cases in which the instantaneous 3-spaces are parallel Euclidean space-like 
hyper-planes not equally spaced due to a linear acceleration. 

A) Rigid non-inertial reference frames with translational acceleration. An 
example are the following embeddings 
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z»(T,a u ) = xZ + e?f(T)+e?a r , 

. 9rT (r,a")=e(^) 2 , g Tr ( T ,a u )=0, g rs ( T , a") = -e 8 rs . 

(3) 

This is a foliation with parallel hyper-planes with normal l^ = e'f = const. 
and with the time- like observer x M (r) = x^ + e% /(r) as origin of the 
3-coordinates. The hyper-planes have translational acceleration x^(t) = 
e^/(r), so that they are not uniformly distributed like in the inertial case 
f(r) = r. 

B) Differentially rotating non-inertial frames without the coordinate sin- 
gularity of the rotating disk. The embedding defining this frames is 

z"(t, a u ) = x»{t) + e» R r s ( T , a) a s x»(t) + e» a r , 

R r s (r,a) = R r s (a t (r,a)) = R r s (F(a) a^r)), 

< F(a) < — !— , dF(a) ^ Q [Moller conditions), 
A a da 



z?(t, a u ) = i"(r) - e{f R r s (r, a) S sw e wuv a u QV ^ a \ 

c 

z?(t, a") = el R k v (r, a) + ^ r) Jr, a) a") , (4) 

where a = \cr\ and R r s {a i (T 1 a)) is a rotation matrix satisfying the asymp- 
totic conditions R r s {i~, a) — ><t_>.oo<^! 9a R t s(t 7 a) -^ a ^oo 0, whose Euler an- 
gles have the expression ai(r, a) = F(a) a,(r), i = 1,2,3. The unit normal 

is I 1 * = e^ = const, and the lapse function is 1 + n(r, er 11 ) = e (z!f 1^ (r, cr M ) = 
ee^ ^(r) > 0. In Eq.(4) one uses the notations Q^(t, a) = i? _1 (r, a) d r R(t, a) 
and [Rr 1 {T,a)d T R{T,a)Y v = S um e mvr " r( c T|g) , with ^ r (r,cr) = J?(r, a) 
n r (r, er) 10 being the angular velocity. The angular velocity vanishes at spa- 
tial infinity and has an upper bound proportional to the minimum of the 
linear velocity vi(t) — x^ l^ orthogonal to the space- like hyper-planes. When 



10 n r (r, <r) defines the instantaneous rotation axis and < Q{t,o) < 
2max I &(t), P(t), j{t) 
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the rotation axis is fixed and J?(r, a) = to = const., a simple choice for the 
function F{a) is F(a) = — 11 • 

1 + " e 2 T 

To evaluate the non-relativistic limit for c — »■ oo, where r = ct with 
i the absolute Newtonian time, one chooses the gauge function F(a) = 
— j 2 a 2 ^c^oo 1 — iiL -r — h 0(c~ 4 ). This implies that the corrections to 

rigidly-rotating non-inertial frames coming from M0ller conditions are of or- 
der 0(c~ 2 ) and become important at the distance from the rotation axis 
where the horizon problem for rigid rotations appears. 

As shown in the first paper in Refs.[21], global rigid rotations are forbid- 
den in relativistic theories, because, if one uses the embedding z^(t,o- u ) = 
x^(t) + e£f R t s (t) o- s describing a global rigid rotation with angular veloc- 
ity fT = £2 t (t), then the resulting g TT {r,o u ) violates M0ller conditions, be- 



cause it vanishes at a = ctr = y x 2 {t) + [ar M (r) £r R r s(T~) (a x i7(r)) r ] 2 

-i„(T)e{fiT a (T)(ff x r)(r)) r j ( a u = aa u , Q r = Q Q r , a 2 = Q 2 = 1). At 

this distance from the rotation axis the tangential rotational velocity becomes 
equal to the velocity of light. This is the horizon problem of the rotating disk 
(the horizon is often named the light cylinder). Let us remark that even if 
in the existing theory of rotating relativistic stars [54] one uses differential 
rotations, notwithstanding that in the study of the magnetosphere of pulsars 
often the notion of light cylinder is still used. 

The search of admissible 3+1 splittings with non-Euclidean 3-spaces is 
much more difficult. The simplest case is the following parametrization of 
the embeddings (1) in terms of Lorentz matrices A A b(t,ct) -^a^oo $b 12 
with A A b(t,0) finite. The Lorentz matrix is written in the form A = BIZ 
as the product of a boost B(t, a) and a rotation 1Z(t, a) like the one in 
Eq.(4) (TT T = 1, K T r = 0, K r s = R r s ). The components of the boost 

are B%{t,<j) = 7 (r,a) = 1/^1 - /3 2 (r, a), B\(r,a) = 7 (t, a) p r (r, a), 

U r s (T,a) = S r s + 2£&( T ,a), with (3 r (T,a) - G(a)0 r (r), whe re (3 r (r) is 

defined by the 4- velocity of the observer (r) = t\ (3 A (r) / \J 1 — /3 2 (t), 

/3 a (t) = (1; /3 r (r)). The M0ller conditions are restrictions on G(cr) -^a^oo 
with G(0) finite, whose explicit form is still under investigation. 

See the second paper of Ref.[21] for the description of the electro-magnetic 
field and of phenomena like the Sagnac effect and the Faraday rotation in this 



11 Nearly rigid rotating systems, like a rotating disk of radius cr , can be described by 
using a function F(a) approximating the step function 8(<r — cr ). 

12 It corresponds to the locality hypothesis of Ref. [20], according to which at each instant 
of time the detectors of an accelerated observer give the same indications as the detectors 
of the instantaneously comoving inertial observer. 
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framework for non-inertial frames. Moreover the embedding (4) has been used 
in the first paper of Ref.[34] on quantum mechanics in non-inertial frames. 



3.3 Congruences of Timelike Observers Associated 
with a 3+1 Splitting 

Each admissible 3+f splitting of space-time allows one to define two associ- 
ated congruences of time-like observers. 

i) The congruence of the Eulerian observers with the unit normal 1^(t, a r ) = 
z^(t, a r ) I a {t 1 er r ) to the 3-spaces embedded in Minkowski space-time as unit 
4-velocity. The world-lines of these observers are the integral curves of the 
unit normal and in general are not geodesies. In adapted radar 4-coordinates 
the contro- variant orthonormal tetrads carried by the Eulerian observers are 

o_A 

l A (r,a r ), A E {a) { Tl cj r ) = {0; 3 e? a) (T,* r )), where 3 ef a) (r,a r ) (a = 1,2,3) are 
triads on the 3-space. 

If 4 V is the covariant derivative associated with the 4- metric 4 Qab{t, <J r ) 
induced by a 3+1 splitting, the equation 

4 V A £h = tlA 3 a B + <7AB+2 3 h AB -uj AB , ( 3 h A B = 4 9ab~£IaIb), (5) 

defines the acceleration a A ( 3 a A I A = 0), the expansion 0, the shear oab = 
o~ba (o~abI B — 0) and the vorticity or twist luab = — <-<Jba (^abI B = 0) 
of the Eulerian observers with ujab = since they are surface-forming by 
construction. They will be useful in GR as shown in Section 7. 

ii) The skew congruence with unit 4-velocity v^(t, cr r ) — z^(t, a r ) v a (t, a r ) 
(in general it is not surface-forming, i.e. it has a non-vanishing vorticity, like 
the one of a rotating disk). The observers of the skew congruence have the 
world-lines (integral curves of the 4-vclocity) defined by a r = const, for every 
t, because the unit 4-velocity tangent to the flux lines x£ o (t) = z^(t, a„) is 

v cr ( T ) = z t ( T > °o ) / \A 4 9tt (t, a r Q ) (there is no horizon problem because it is 
everywhere time- like in admissible 3+1 splittings). They carry contro- variant 
orthonormal tetrads, given in Ref.[41], not adapted to the foliation, connected 
in each point by a Lorentz transformation to the ones of the Eulerian observer 
present in this point. 
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3-4 Parametrized Minkowski Theories 

In the global non-inertial frames of Minkowski space-time it is possible to 
describe isolated systems (particles, strings, fields, fluids) admitting a La- 
grangian formulation by means of parametrized Minkowski theories [23, 24], 
[21]. 

The existence of a Lagrangian, which can be coupled to an external 
gravitational held, makes possible the determination of the matter energy- 
momentum tensor and of the ten conserved Poincare generators P M and J^ v 
(assumed finite) of every configuration of the isolated system. 

First of all one must replace the matter variables of the isolated system 
with new ones knowing the clock synchronization convention defining the 
3-spaccs S T . For instance a Klein-Gordon field <fr(x) will be replaced with 
4>{t, <J r ) = <P(z(t, <y r ))\ the same for every other field. Instead for a rela- 
tivistic particle with world-line x M (r) one must make a choice of its energy 
sign: then the positive- (or negative-) energy particle will be described by 
3-coordinates rf(r) defined by the intersection of its world- line with E T : 
x p (t) = z^(t, J7 r (r)). Differently from all the previous approaches to relativis- 
tic mechanics, the dynamical configuration variables are the 3-coordinatcs 
77 r (r) and not the world- lines x m (t) (to rebuild them in an arbitrary frame 
one needs the embedding defining that frame). This fact eliminates the pos- 
sibility to have time-like excitations in the spectrum of rclativistic bound 
states: inside each 3-space only space-like correlations among the particles 
are possible. 

Then one replaces the external gravitational 4-metric in the coupled La- 
grangian with the 4-metric 5ab(t, er r ), which is a functional of the embed- 
ding defining an admissible 3+1 splitting of Minkowski space-time, and the 
matter fields with the new ones knowing the instantaneous 3-spaces S T . 

Parametrized Minkowski theories are defined by the resulting Lagrangian 
depending on the given matter and on the embedding z^(T,a r ). The re- 
sulting action is invariant under tto take into account of fermion fieldshc 
frame-preserving diffeomorphisms r h- > r (r, cr u ), o~ r i-> a r (cr u ) firstly intro- 
duced in Ref.[25]. As a consequence, there are four first-class constraints with 
exactly vanishing Poisson brackets (an Abelianized analogue of the super- 
Hamiltonian and super-momentum constraints of canonical gravity) deter- 
mining the momenta conjugated to the embeddings in terms of the matter 
energy-momentum tensor. This implies that the embeddings z^{t,o~ t ) arc 
gauge variables, so that all the admissible non-inertial or inertial frames are 
gauge equivalent, namely physics does not depend on the clock synchroniza- 
tion convention and on the choice of the 3-coordinates cr r : only the appear- 
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anccs of phenomena change by changing the notion of instantaneous 3-space 

13 

Even if the gauge group is formed by the frame-preserving diffcomor- 
phisms, the matter energy-momentum tensor allows the determination of 
the ten conserved Poincare generators P M and J^ v (assumed finite) of ev- 
ery configuration of the system (in non-incrtial frames they are asymptotic 
generators at spatial infinity like the ADM ones in GR). 

As an example one may consider N free scalar particles with masses rrii 
and sign of the energy 77, = d=, whose world- lines are identified by the conhg- 
urational variables ?7[(t): x?(t) — z^(t,7]1(t)), i — 1,..,N. In parametrized 
Minkowski theories they are described by the following action depending on 
the configurational variables f?f (t) and z^(t, a r ) 



S = J dTd 3 <r£(T 7 <7 u )= J drL(r), 



N 



£(r,a") = -^5 3 K-^(r)) 



1=1 



mem yje [ 4 3tt (t, <t«) |2V(t, <t«) r,l(T) + 4 ffrs (r, <t«) ^(t) 77? (r)]. 

(6) 

The resulting canonical momenta Ki r (r) — d - , Ph{t, <J u ) = — e q ^{rxr^)) 
satisfy the Poisson brackets {^[(t), Kj s (r)} = —6% {z^(t, a w ), p v {t, o ")} = 
— <5£f <5 3 (cr" — <7 u ). The Poincare generators and the energy-momentum tensor 
of this system are (h rs — —e-f rs with "f ru4 9us = 5 r s ] 7 = —edet 4 g rs ) 

P» = j d 3 a(f{T,a u ), = J d z a{z^p v - z"pV)( T ,a u ), 

T AB (t rr u ) — - - - mfiv _ u v rpAB 

[T ' a) - y/-det*gc D (T,o«) 6*gAB(T t o~y ~ Z * Zb1 ' 



\/i{t, <? u ) 



Vi \J m i c2 + h rs (r, <7 u )n ir (T) k is (t), 



13 In the first paper of Ref.[34] there is the definition of parametrized Galilei theories, non 
rclativistic limit of the parametrized Minkowski theories. Also the inertial and non-incrtial 
frames in Galilei space-time are gauge equivalent in this formulation. 
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N 



T ±r (r,a u ) 



T rs (T,<T U ) 




(7) 



^/mfc 2 +h vw (T, <T u )K iv (T) K iw (T)' 



The four first-class constraints implying the gauge nature of the embedding 
and the gauge equivalence of the description in different non-inertial frames 
are 



The same description can be given for the Klein-Gordon [55] and Dirac 
[56] fields and for the electro-magnetic field [21]. 

To describe the physics in a given admissible non-inertial frame described 
by an embedding Zp(r, a u ) one must add the gauge- fixings z^(t 7 <t u ) — 



3. 5 The Instant Form of Dynamics in the Inertial 
Rest Frames and the Problem of the Relativistic 
Center of Mass 

If one restricts himself to inertial frames, one can define the inertial rest-frame 
instant form of dynamics for isolated systems by choosing the 3+1 split- 
ting corresponding to the intrinsic inertial rest frame of the isolated system 
centered on an inertial observer: the instantaneous 3-spaces, named Wigner 
3- spaces due to the fact that the 3-vectors inside them are Wigner spin- 1 3- 
vectors [23, 24], are orthogonal to the conserved 4-momcntum (assumed 
time-like, e P 2 > 0) of the configuration. 

In this framework one can give the final solution to the old problem of the 
relativistic extension of the Newtonian center of mass of an isolated system. 
In its rest frame there are only three notions of collective variables, which 
can be built by using only the Poincare generators : 

the canonical non-covariant Newton- Wigner center of mass (or center of 
spin) x^(t), 

the non-canonical covariant Fokker-Pryce center of inertia y(r) 
the non-canonical non-covariant M0ller center of energy i? Ai (r). 
While Y^(t) is a 4- vector, x M (r) and i? M (r) are not 4- vectors. All of them 
tend to the Newtonian center of mass in the non- relativistic limit. Since the 



p„(T,<r u )- V7(t,<7«) l^T ±± -z r ^h rs T ±s (t,O»0. 



(8) 
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Poincare generators know the whole E T , they and therefore also these three 
collective variables are non-local quantities: as a consequence they are non 
measurable with local means [21, 26, 28, 33]. 

If one centers the inertial rest frame on the world- line of the Fokkcr-Planck 
center of inertia thought as an inertial observer, then the corresponding em- 
bedding has the expression [23, 24, 28] 

z» v (T,<T) = Y»(T)+e?(h)o- r , (9) 
where Y^(t) is the Fokkcr-Prycc ccntcr-of-incrtia 4- vector, h — P/Ve P 2 

o 

and t 11 a=v (h) = L^A=iy(P, P) are the columns of the standard Wigner boost 

o ft 

for time-like orbits sending P» = VTP 2 (VI + h 2 : h) to P = VeP 2 (1; 0). 
Their expression is e£(h) = = + h 2 ; h) and ef{h) = (h r ;8 l r + 

— — ) as shown in Appendix B of Ref.[57l. 
i+VT+hV 

As shown in Ref.[26, 28, 29, 33], the three collective variables can be ex- 
pressed as known functions of the Lorentz-scalar rest time r = cT s = h ■ x = 
h ■ Y = h ■ R, of canonically conjugate Jacobi data (frozen Cauchy data) 
z = Mcx NW (0) and h = P/Mc 14 , of the invariant mass Mc = V 'eP 2 of 
the system and of its rest spin S. 

While the world-line of the non-canonical covariant external Fokker-Pryce 
4-center of inertia is 



Y»(t)=z&{t,0)=(x°(t);Y(t)) = 

/ r -o hz. z h z. 1 Sxh \ 

= VI + h 2 r + -rr- ; — + (t + — - h + , 

V Mc h Mc Mc' Mc(l + VT+^) J 

(10) 

the pseudo- world-line of the canonical non-covariant external 4-center of mass 
is (cr = from Ref.[26]) 

i"(r) = (V(r);*(r)) = ^(r,*) = r"(r) + (0; " S>< i^ ) = 

v y v Mc(l + Vl + h r 



14 Their Poisson brackets a,re{z' t , hi} = & l i . xj^yv{r) is the standard Newton- Wigner non- 
covariant 3-position, classical counterpart of the corresponding position operator; the use of 
z avoids to take into account the mass spectrum of the isolated system in the description of 
the center of mass. The non-covariance of z under Poincare transformations (a, A) has the 

following form [33, 57] z i ^ z i = (a^ - ^ ^ X° j j z' + ( A\ - ) C^ 1 a) M - 



21 



Luca Lusanna 



(ii) 



The world-lines of the positive-energy particles are parametrized by the 
Wigner 3-vectors tj^t), i = 1,2, .., N, and are given by 



The world-lines (r) of the particles are derived (interaction-dependent) 
quantities. Also the standard particle 4-momenta arc derived quantities, 
whose expression is pf (r) = e^(h) Kf(r) = W ^mfc 2 + k 2 (t) — e£(h) K, r (r) 
with epf = mfc 2 in the free case. 

In the case of interacting particles the reconstruction of the world-lines 
requires a complex interaction-dependent procedure delineated in Ref. [29], 
where there is also a comparison of the present approach with the other for- 
mulations of rclativistic mechanics developed for the study of the problem of 
relativistic bound states. See Ref. [21] for the extension to non-inertial frames. 

In general the world-lines xf (r) do not satisfy vanishing Poisson brackets 
(they are relativistic predictive coordinates, see Ref. [29]): already at the clas- 
sical level a non- commutative structure emerges due to the Lorentz signature 
of the space-time [33] . 

In each Lorentz frame one has different pseudo-world-lines describing 
and x^: the canonical 4-center of mass x^ lies in between and in 
every (non rest)-frame. As discussed in Subsection IIF of Ref. [28] , this leads 
to the existence of the Miller non-covariance world-tube, around the world- 
line F M of the covariant non-canonical Fokker-Pryce 4-center of inertia Y^ . 
The invariant radius of the tube is p = y/—e W 2 /p 2 = |S|/Ve P 2 where 
(W 2 = —P 2 S 2 is the Pauli-Lubanski invariant when e P 2 > 0). This classical 
intrinsic radius is a non-local effect of Lorentz signature absent in Euclidean 
spaces and delimits the non-covariance effects (the pseudo- world-lines) of the 
canonical 4-center of mass x^ 15 . They are not detectable because the M0llcr 
radius is of the order of the Compton wave-length: an attempt to test its 
interior would mean to enter in the quantum regime of pair production. The 
M0ller radius p is also a remnant of the energy conditions of general relativity 
in flat Minkowski space-time [23] . 

Finally Eqs.(7) can be used to extend the multipolar expansions of Ref. [58] 
to this framework for relativistic isolated systems as it is shown in the third 
paper of Refs.[26]. 



15 In the rest-frame the world-tube is a cylinder: in each instantaneous 3-space there is 
a disk of possible positions of the canonical 3-center of mass orthogonal to the spin. In 
the non-relativistic limit the radius p of the disk tends to zero and one recovers the non- 
rclativistic center of mass. 



x?(t) = z^(r, Vi (r)) = Y»(T)+e?(T)r)l(T). 



(12) 
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3.6 The Description of Isolated Systems in the Rest 
Frame and their Poincare Generators 

In the inertial rest frame of an isolated system Eqs.(7) are the starting point 
to get the explicit form of its Poincare generators, in particular of the Lorentz 
boosts, which, differently from the Galilei ones, are interaction dependent. 

As shown in Ref.[28], every isolated system (i.e. a closed universe) can be 
visualized as a decoupled non-covariant collective (non-local) pseudo-particle 
(the external center of mass), described by the frozen Jacobi data z, h, car- 
rying a pole-dipole structure, namely the invariant mass M c (the monopole) 
and the rest spin S (the dipole) of the system, and with an associated external 
realization of the Poincare group 16 



J ij = z l h? - z j ti + e ijk S k , K l = J' 



r T2 i {Sxhy 

= -Vl + hz l + 



1 + Vl + h 2 ' 
(13) 



The universal breaking of Lorentz covariance is connected to this decou- 
pled non-local collective variable and is irrelevant because all the dynam- 
ics of the isolated system leaves inside the Wigner 3-spaces and is Wigner- 
covariant. The invariant mass and the rest spin are built in terms of the 
Wigner-covariant variables of the given isolated system (tj^t) and k^(t)) 
inside the Wigner 3-spaces [21, 26, 28, 33]. 

In each Wigner 3-space E T there is a unfaithful internal realization of the 
Poincare algebra, whose generators are built by using the energy-momentum 
tensor (7) of the isolated system. While the internal energy and angular 
momentum are Mc and S respectively, the internal 3-momentum vanishes: 
it is the rest frame condition. Also the internal (interaction dependent) 
Lorentz boost vanishes: this condition identifies the covariant non-canonical 
Fokker-Pryce center of inertia as the natural inertial observer origin of the 
3-coordinates a r in each Wigner 3-space. 

For N free particles the internal Poincare generators have the following 
expression 



16 The last term in the Lorentz boosts induces the Wigner rotation of the 3-vectors inside 
the Wigner 3-spaces. 
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Mc=- £ {mt) = \J m\ c 2 + k 
c ;=i 

N 

'P(int) — J]] K i ~ 0, 

i=l 

AT 

S = (int) — J]] x K ii 



i=l 

N 



JC {int) = - c 2 + k 2 « 0. (14) 



i=i 



Since one is in an instant form of the dynamics, in the interacting case 
only Mc and /C( int ) become interaction dependent. 

The three pairs of second-class (interaction dependent) constraints T(i n t) ~ 
0, 7C( int ) w 0, eliminate the internal 3-center of mass and its conjugate mo- 
mentum inside the Wigner 3-spaces 17 : this avoids a double counting of the 
collective variables (external and internal center of mass). As a consequence 
the dynamics inside the Wigner 3-spaces is described in terms of internal 
Wigner-covariant relative variables. In the case of N relativistic particles one 
defines the following canonical transformation [33] (see Ref. [26] for other vari- 
ants) (m = J2iLi m i) 



N N 

— Vi, K+=V (int )=2^Ki, 
i=l i=l 
N N 
Pa = V^V ^ laiVi, ^a = -7=j^2r ai Ki, O = 1, .., JV - 1, 

i=l * i=l 



N-l N-l 



Vi = 



V+ + ~7= Fal P *> Kl = ~ K + + ^ Y Tai^a, ( 15 ) 

ViV o-l m a=l 

with the following canonicity conditions 18 



17 One can show [26, 28] that one has 7C(j nt \ = — M R+, where R+ is the internal M0ller 
3-center of energy inside the Wigner 3-spaces. The rest frame condition "Pn n £\ ~ implies 

RJ q + ?i V+> where q + is the internal 3-center of mass and y + the internal Fokker- 
Pryce 3-center of inertia. 

18 Eqs.(15) describe a family of canonical transformations, because the 7 a ;'s depend on 
5 (AT — 1)(N — 2) free independent parameters. 
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^ ' lai — 0, ^ ' lai Ibi 3ab, ^ ' lai laj $ij ,y i 

i=l i=l a=l 

N m ■ 1 N 

Fax lai / J lak ; lai ^ai "77 / , ^ak ; 

N N N-l 

/ , 1 at — u ; / , fat 1 bi — u abj / , Jai 1 ai — u ij 

^ m f—' ^— ; m 

i— 1 2—1 a=l 



(16) 



Since Eqs.(14) imply K+(r) = P (i „ t) w and r/ + (r) w /_ t (p (r), 7T a (r)) 
due to /C( int ) ~ 0, the invariant mass Mc and the rest spin S become func- 
tions only of the N-l pairs of relative canonical variables. 

As a consequence, Eqs.(10), (12) and (15) imply that the world-lines (r) 
can be expressed in terms of the Jacobi data z, h, and of the relative vari- 
ables P q (t), 7r a (r), a = 1,..,N — 1. See Ref.[55] for the collective and rel- 
ative variables of the Klein-Gordon field and the second paper in Ref.[28] 
for such variables for the electro-magnetic field in the radiation gauge. For 
these systems one can give for the first time the explicit closed form of the 
interaction-dependent Lorentz boosts. 

One finds that disregarding the unobservable external center of mass all 
the dynamics is described only by relative variables: this is a form of weak 
relationism without the heavy foundational problem of approaches like the 
one in Ref.[59]. 



The non-relativistic limit of this description [33] is Newton mechanics with 
the Newton center of mass decoupled from the relative variables and moreover 
after a canonical transformation to the frozen Hamilton- Jacobi description 
of the center of mass. 



An important remark is that the internal space of relative variables is 
independent from the reference inertial frame used for the description of 
the isolated system. As shown in Ref.[33], the formalism is built in such a 
way that the Wigner rotation induced on the relative variables by a Lorentz 
transformation connecting two reference inertial frames is the identity, i.e. 
the space of the relative variables in an abstract internal space insensitive to 
Lorentz transformations carried by the external center of mass (or in a more 
covariant description carried by the Fokker-Pryce center-of-inertia 4-vector 
origin of the embedding (9)). 

Finally in Ref. [21] there is the extension of the formalism to admissible 
non-inertial rest frames, where P M is orthogonal to the asymptotic space-like 
hyper-planes to which the instantaneous 3-spaces tend at spatial infinity. In 
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these non-inertial rest frames the internal Poincare generators are asymptotic 
(constant of the motion) symmetry generators like the asymptotic ADM ones 
in the asymptotically Minkowskian space-times. 
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4 Implications for Relativistic Mechanics and Classical 
Field Theory in Special Relativity and the 
Multi- Temporal Quantization Approach 

In the rest-frame instant form of the dynamics it has been possible to find 
the explicit form of the internal Poincare generators (in particular of the 
interaction-dependent invariant mass and Lorentz boosts) not only for the 
Klein-Gordon [55] and Dirac [31, 56] fields, but also for the electro-magnetic 
field in the radiation gauge (the only one suitable for the Shanmugadhasan 
canonical transformations of constraint theory [60]) [21, 28], for relativistic 
fluids [41, 42], spinning particles [27, 31] and for massless particles, the Nambu 
string and the two- level atom [32]. 

In this Section some other developments in SR will be reviewed. 



4-1 Relativistic Atomic Physics 

Standard atomic physics [61] is a semi-relativistic treatment of quantum 
electro-dynamics (QED) in which the matter fields are approximated by 
scalar (or spinning) particles, the relevant energies are below the threshold of 
pair production and the electro-magnetic field is described in the Coulomb 
gauge at the order 1/c. 

In Rcfs.[27, 28] a fully relativistic formulation of classical atomic physics in 
the rest-frame instant form of dynamics was given with the electro-magnetic 
field in the radiation gauge and with the electric charges Qi of the positive- 
energy particles being Grassmann-valucd (Qf = 0, Qi Qj — Qj Qi for i ^ j) 
to regularize the electro-magnetic self-energies on the world-lines of parti- 
cles. In the language of QED this is both a ultraviolet regularization (no loop 
contributions) and an infrared one (no brehmstrahlung) , so that only the one- 
photon exchange diagram contributes and its static and non-static effects are 
replaced by potentials in a formulation based on the Cauchy problem. There- 
fore the starting point is a parametrized Minkowski theory with N charged 
positive-energy particles mutually interacting with a Coulomb potential and 
coupled to a dynamical transverse electro-magnetic field described by the 
canonical variables A±(r,a r ) and tv±(t, a r ) = E±(r,a r ). 

In the first paper of Rcf . [27] (the second paper is devoted to spinning par- 
ticles) it is shown that the use of the Lienard-Wiechert solution (see the third 
paper in Ref.[27]) with "no incoming radiation field" allows one to arrive at 
a description of N charged particles dressed with a Coulomb cloud and mu- 
tually interacting through the Coulomb potential augmented with the full 
relativistic Darwin potential. This happens independently from the choice of 
the Green function (retarded, advanced, symmetric,..) due to the Grassmann 
regularization. The quantization allows one to recover the standard instan- 
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taneous approximation for relativistic bound states, which till now had only 
been obtained starting from QED (either in the instantaneous approxima- 
tions of the Bethc-Salpeter equation or in the quasi-potential approach). In 
the case of spinning particle the relativistic Salpeter potential was identified. 

Moreover in Ref. [28] it is shown that by using the previous results one 
can find a canonical transformation from the canonical basis »7j(r), k»(t), 
A±(t, cr r ), 7t^(t, CT r ), in which the internal Poincare generators have the ex- 
pression in the case N=2 (B = d x A±, c(<r) = — l/An \a\) 



r£ 4?r I Vi(r)-rij(r) I 2 J 

N 1 

P(int) =^2*i(T) + - / d 3 a [ir± x B](t, <j) m 0, 
i=i c J 

5 r = E (r)i(T) x Ki (T)) r + - I d 3 a(<r * (W xB])V,ff), 

i—1 



^»*) = -E^( T )v m?c2 + ( K4(T) "T A±(T ' ^ (t)) ) + 

AT 1..JV _ 



i ^ r — ^ ^ r j ^ 

-^ r ( T ) c( r? 4 (r)-^.(T))] + 

+ Q,y d 3 <7 71"^ (r, <t) c(er — Tfc(r))] - ^ (tiI + £ 2 )(r, W,) 

to a new canonical basis r)i(r), Ki(r), Ai ra( i(r,(T r ), 7rj_ ra d(r, a r ) so that 
in the rest frame there is a decoupled free radiation transverse field and a 
system of charged particles mutually interacting with Coulomb plus Darwin 
potential. See the first paper in Ref. [27] for the explicit form of the relativistic 
Darwin potential. The new internal Poincare generators in the N=2 case are 



From Clock Synchronization 31 

+ 2 J d3(J ( n ±rad + B lad) ( T ^ (T ) = ^matter + £rad, 
2 j - 

P(int) = ^ «z(t) + ~ / d'V ( 7T irad X B rad J (T, Cr) = T matter + "Prad ~ 0, 
i=l CJ 

S = ^ f) l X ki + - J d 3 (T <J X ^7T ira(J X B rad ^j ("T, Cr) = S 'matter + S ra d, 



2 



K(int) = - ^ »7 4 y m? c 2 + k - - 



i=i 



- / 1 9 /Ci2(ki,k 2 ,Pi 2 ) o 4 ( C ^ 

IQiQ2\, Kl "U 3Pi2 -2A ±S2 (K 2 ,p 12; 

*7i 



2 c 



2 C 2 + K 2 



+ *?2 " 



2 9 i ' 2 



<3lQ2 / / 2 2 , "2 9 / 2 2,-2 O U r ~ - x 
" [\J m l c + K 1 + V m 2 C + K2 d^) ^-12( K l> K 2,Pl 2 ) - 

?1 Q2 /"^ / TT±Sl(o- - 7T±g 2 (<T ~ r? 2 ^2) \ 

4?rc J V |cr — T> 2 j \<r-fli\ ' 



QiQ 



c 



J d 3 aa (tv ±S i(<t - Vi, «i) • 7r_L,s2 (er -f] 2 ,k 2 ) + 
+ B S l{(T - Th,Ki) • B S 2(cr - ^2^2)) - 

- Y~ c j d3<J i^^rad + B lad) ( T ; °") = ^-matter + K-rad ~ 0. (18) 

The only restriction on the two decoupled systems is the elimination of their 
overall internal 3-center of mass inside the Wigner 3-spaces. Therefore, at 
the classical level there is a way out from the Haag theorem forbidding the 
existence of the interaction picture in QED, so that there is no unitary evo- 
lution based on interpolating fields from the "in" states to the "out" ones in 
scattering processes. While the extension of these results to the non-inertial 
rest frame is done in Rcf. [21], the quantization of this framework is under 
investigation. 

In the first paper of Ref.[32] there is the formulation in the rest- frame 
instant form of the relativistic quark model in the radiation gauge for the 
SU(3) Yang-Mills fields with scalar quarks having Grassmann-valued color 
charges. While in Eq.(lOl) of that paper there is the rest-frame condition, 
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in Eqs.(97) there is the invariant mass Mc 2 for a quark- ant iquark system. 
In it the electro- magnetic Coulomb potential of Eq.(17) is replaced with a 
potential, given in Eq.(95), depending on the color transverse vector po- 
tential through the Green function of the SU(3) covariant divergence. The 
non-linearity of the problem does not allow to evaluate a Lienard-Wiechert 
solution and to find the analogue of Eqs.(18) . 



4-2 Relativistic Kinetic Theory and Relativistic 
Micro- Canonical Ensemble 

In the rest-frame instant form of dynamics it is also possible to give a finally 
consistent treatment of relativistic kinetic theory and relativistic statistical 
mechanics [30] . In particular one can give a definition of the relativistic micro- 
canonical ensemble for an isolated system of N interacting particles with fixed 
internal energy £ and rest spin S only in terms of the internal Poincare 
generators in the Wigner 3-spaces by means of the partition function (V is 
the volume defined by the function x(V) vanishing outside it) 

1..JV 1..JV 

Z(£,S,V,N) = — / J] d\ x (V) I J] d 3 K,S(Mc 2 -£) 

— , - 

Also it extension to non-incrtial rest frames can be given by using the 
results of Ref.[21] with the result that notwithstanding the presence of long- 
range inertial forces one has still an equilibrium distribution. 



4-3 Relativistic Quantum Mechanics and Relativistic 
Entanglement 

A new formulation of relativistic quantum mechanics in the Wigner 3-spaces 
of the inertial rest frame is developed in Ref. [33] in absence of the electro- 
magnetic field. It englobes all the known results about relativistic bound 
states (absence of relative times) and avoids the causality problems of the 
Hegerfeldt theorem [62] (the instantaneous spreading of wave packets). 

In it one quantizes the frozen Jacobi data z and h of the canonical non- 
covariant decoupled external center of mass and the relative variables in the 
Wigner 3-spaces. Since the center of mass is decoupled, its non-covariance is 
irrelevant: like for the wave function of the universe, who will observe it? 
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The resulting Hilbert space has the following tensor product structure: 
HJ <g> H re i, where H conitHJ is the Hilbert space of the exter- 
nal center of mass (in the Hamilton- Jacobi formulation due to the use of 
frozen Jacobi data) while H re i is the Hilbert space of the relative variables in 
the abstract internal space living in the Wigner 3-spaces. While at the non- 
relativistic level this presentation of the Hilbert space is unitarily equivalent 
to the tensor product of the Hilbert spaces Hi of the individual particles 
H = H\ <g) H2 ® this is not true at the relativistic level. 

If one considers two scalar quantum particles with Klein-Gordon wave 
functions belonging to Hilbert spaces "H x ° , in the tensor-product Hilbert space 
(Hi) x ° ® (^2)1° ® ••• there is no correlation among the times of the particles 
(their clocks are not synchronized) so that in most of the states there are 
some particles in the absolute future of the others. As a consequence the two 
types of Hilbert spaces lead to unitarily inequivalent descriptions and have 
different scalar products (compare Refs. [33] and [57]). 

As a consequence, at the relativistic level the zeroth postulate of non- 
relativistic quantum mechanics does not hold: the Hilbert space of composite 
systems is not the tensor product of the Hilbert spaces of the sub-systems. 
Contrary to Einstein's notion of separability (separate objects have their 
independent real states) [63] one gets a kinematical spatial non-separability 
induced by the need of clock synchronization for eliminating the relative times 
and to be able to formulate a well-posed relativistic Cauchy problem. 

Moreover one has the non-locality of the non-covariant external center of 
mass which implies its non-measurability with local instruments 19 . While its 
conjugate momentum operator must be well defined and self-adjoint, because 
its eigenvalues describe the possible values for the total momentum of the 
isolated system (the momentum basis is therefore a preferred basis in the 
Hilbert space), it is not clear whether it is meaningful to define center-of- 
mass wave packets. 

These non-locality and kinematical spatial non- separability are due to the 
Lorentz signature of Minkowski space-time and this fact reduce the relevance 
of quantum non-locality in the study of the foundational problems of quantum 
mechanics [63] which have to be rephrased in terms of relative variables. 

The quantization defined in Rcf. [33] leads to a first formulation of a the- 
ory for relativistic entanglement, which is deeply different from the non- 
relativistic entanglement due to these kinematical non-locality and spatial 
non-separability. To have control on the Poincare group one needs an iso- 
lated systems containing all the relevant entities (for instance both Alice and 
Bob) of the experiment under investigation and also the environment when 
needed. One has to learn to reason in terms of relative variables adapted to 



19 In Ref. [64] it was shown that the quantum Newton- Wigner position should not be a 
self-adjoint operator, but only a symmetric one, with an implication of bad localization. 
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the experiment like molecular physicists do when they look to the best sys- 
tem of Jacobi coordinates adapted to the main chemical bonds in the given 
molecule. This theory has still to be developed together with its extension to 
non-inertial rest frames. 



4-4 Multitemporal Quantization in Non-inertial 
Frames 

This quantization of relativistic mechanics can be extended to the class of 
global non-inertial frames with space-like hyper-planes as 3-spaces and differ- 
entially rotating 3-coordinates defined in Ref.[21] by using the multi-temporal 
quantization approach developed in Rcf. [65] . 

As shown in Rcf. [34] , in this type of quantization one quantizes only the 3- 
coordinates r\\ (r) of the particles and not the inertial effects (like the Coriolis 
and centrifugal ones): they remain c-numbers describing the appearances of 
phenomena. The known results in atomic and nuclear physics are reproduced. 



4-5 Open Problem 

The main open problem in SR is the quantization of fields in non-inertial 
frames due to the no-go theorem of Ref.[66] showing the existence of ob- 
structions to the unitary evolution of a massive quantum Klein-Gordon field 
between two space-like surfaces of Minkowski space-time. It turns out that 
the Bogoljubov transformation connecting the creation and destruction op- 
erators on the two surfaces is not of the Hilbcrt-Schmidt type, i.e. that the 
Tomonaga-Schwinger approach in general is not unitary. One must reformu- 
late the problem using the nice foliations of the admissible 3+1 splittings 
of Minkowski space-time and to try to identify all the 3+1 splittings allow- 
ing unitary evolution. This will be a prerequisite to any attempt to quantize 
canonical gravity taking into account the equivalence principle (global iner- 
tial frames do not exist) with the further problem that in general the Fourier 
transform does not exist in Einstein space-times. 
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5 Non-Inertial Frames in Asymptotically Minkowskian 
Einstein Space-Times and ADM Tetrad Gravity 

After this description of SR induced by the metrology-oriented problem of 
clock synchronization, one has to face the same problems in the globally hy- 
perbolic, topologically trivial, asymptotically Minkowskian space-times with- 
out super-translations 20 of GR. As shown in the first paper of Ref.[17], in 
the chosen class of space-times the ten strong asymptotic ADM Poincare gen- 
erators P^dm > J adm (they are fluxes through a 2-surface at spatial infinity) 
are well defined functionals of the 4-metric fixed by the boundary conditions 
at spatial infinity. 

While in SR Minkowski space-time is an absolute notion, in Einstein GR 
also the space-time is a dynamical object [35] and the gravitational field 
is described by the metric structure of the space-time, namely by the ten 
dynamical fields A g IJtv {x) (x^ are world 4-coordinates) . The 4-metric 4 g fl „(x) 
tends in a suitable way to the flat Minkowski 4-metric 4 ?? MI , at spatial infinity 
[17]: having an asymptotic Minkowskian background the usual splitting of the 
4-metric in the bulk in a background plus perturbations in the weak field 
limit can be avoided as shown in Section VII. 

The ten dynamical fields A g^ v {x) are not only a (pre)potential for the 
gravitational field (like the electro-magnetic and Yang-Mills fields are the 
potentials for electro-magnetic and non-Abelian forces) but also determines 
the chrono- geometrical structure of space-time through the line element ds 2 = 
4 g^i, dx^ dx" '. Therefore the 4-metric teaches relativistic causality to the other 
fields: it says to massless particles like photons and gluons which are the 
allowed world-lines in each point of space-time. This basic property is lost 
in every quantum field theory approach to gravity with a fixed background 
4-metric 21 . 

In these space-times one can define global non-inertial frames by using 
the same admissible 3+1 splittings, centered on a time-like observer, and 
the observer-dependent radar 4-coordinates a A — (r;a r ) employed in SR. 
This will allow to separate the inertial (gauge) degrees of freedom of the 
gravitational field (playing the role of inertial potentials) from the dynamical 
tidal ones at the Hamiltonian level. 



These space-times must also be without Killing symmetries, because, otherwise, at the 
Hamiltonian level one should introduce complicated sets of extra Dirac constraints for each 
existing Killing vector. 

21 The ACES mission of ESA [67] will give the first precision measurement of the grav- 
itational red-shift of the gcoid, namely of the 1/c 2 deformation of Minkowski light-cone 
caused by the geo-potential. In every quantum field theory approach to gravity, where the 
definition of the Fock space requires the use of the Fourier transform on a fixed back- 
ground space-time with a fixed light-cone, this is a non-pcrturbative effect requiring the 
re-summation of the perturbative expansion. 
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In GR the dynamical fields are the components A g liv (x) of the 4-metric 
and not the embeddings x^ — z^{T,a r ) dehning the admissible 3+1 split- 
tings of space-time like in the parametrized Minkowski theories of SR. Now 
the gradients z^(t, <r r ) of the embeddings give the transition coefficients 
from radar to world 4-coordinates, so that the components 4 gAB(i~, cr r ) = 
Za(t, o r ) Zg(r, a r ) 4 g l _ Ll/ (z(T, <r r )) of the 4-metric will be the dynamical fields 
in the ADM action. Like in SR the 4- vectors z^(t, cr r ), tangent to the 3-spaces 
E T , are used to define the unit normal i^(r, a r ) — z^(t, a r ) I a (t, a r ) to S T , 
while the 4- vector z£(t, <r r ) has the lapse function as component along the 
unit normal and the shift functions as components along the tangent vectors. 

Since the world-line of the time-like observer can be chosen as the origin 
of a set of the spatial world coordinates, i.e. x M (r) = (x°(t);0), it turns 
out that with this choice the space-like surfaces of constant coordinate time 
x°(t) = const, coincide with the dynamical instantaneous 3-spaces S T with 
t = const.. By using asymptotic flat tetrads = 5% 8 T A + 8? 8\ (with e A 
denoting the inverse flat cotetrads) and by choosing a coordinate world time 
x°(t) = x° + e°r = x° + t, one gets the following preferred embedding 
corresponding to these given world 4-coordinates 

x» = z»(t, o r ) = x»(t) + e£ a r = 8% x° + t\ a A . (20) 
This choice implies z^(t, a r ) — and 4 g^(x — z(r,a r )) — e A 4 gAB(i~, a r ). 

As shown in Ref.[35], the dynamical nature of space-time implies that 
each solution (i.e. an Einstein 4-geometry) of Einstein's equations (or of the 
associated ADM Hamilton equations) dynamically selects a preferred 3+1 
splitting of the space-time, namely in GR the instantaneous 3-spaces are 
dynamically determined in the chosen world coordinate system. Eq.(20) can 
be used to describe this 3+1 splitting and then by means of 4-diffeomorphisms 
the solution can be written in an arbitrary world 4-coordinate system in 
general not adapted to the dynamical 3+1 splitting. This gives rise to the 
4-geometry corresponding to the given solution. 

To define the canonical formalism the Einstcin-Hilbcrt action for metric 
gravity (depending on the second derivative of the metric) must be replaced 
with the ADM action (the two actions differ for a surface tern at spatial in- 
finity). As shown in the first paper of Rcfs.[17], the Legendre transform and 
the definition of a consistent canonical Hamiltonian require the introduction 
of the DeWitt surface term at spatial infinity: the final canonical Hamilto- 
nian turns out to be the strong ADM energy (a flux through a 2-surface at 
spatial infinity), which is equal to the weak ADM energy (expressed as a 
volume integral over the 3-space) plus constraints. Therefore there is not a 
frozen picture but an evolution generated by a Dirac Hamiltonian equal to 
the weak ADM energy plus a linear combination of the first class constraints. 
Also the other strong ADM Poincare generators are replaced by their weakly 
equivalent weak form P^dm^ ^adm- 
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In the first paper of Ref.[17] it is also shown that the boundary conditions 
on the 4-metric required by the absence of super-translations imply that the 
only admissible 3+1 splittings of space-time (i.e. the allowed global non- 
inertial frames) are the non-inertial rest frames: their 3-spaces are asymp- 
totically orthogonal to the weak ADM 4-momentum. Therefore one gets 
P r ADM ~ as the rest-frame condition of the 3-universe with a mass and a 
rest spin fixed by the boundary conditions. Like in SR the 3-universe can be 
visualized as a decoupled non-covariant (non-measurable) external relativis- 
tic center of mass plus an internal non-inertial rest-frame 3-space containing 
only relative variables (see the first paper in Ref.[40]). 



5.1 The Parametrization of Tetrads for ADM Tetrad 
Gravity 

To take into account the coupling of fermions to the gravitational field metric 
gravity has to be replaced with tetrad gravity. This can be achieved by de- 
composing the 4-metric on cotetrad fields (by convention a sum on repeated 
indices is assumed) 

4 9AB(r,a r ) = E^\r,a r )\ am E^(r,a r ), (21) 

by putting this expression into the ADM action and by considering the re- 
sulting action, a functional of the 16 fields (r, a r ), as the action for ADM 
tetrad gravity. In Eq.(21) (a) are flat indices and the cotetrad fields E A °^ are 
the inverse of the tetrad fields E A s , which are connected to the world tetrad 

fields by E» a) {x) = ^(r, cr r ) (z(r, cr r )) by the embedding of Eq.(20). 

This leads to an interpretation of gravity based on a congruence of time- 
like observers endowed with orthonormal tetrads: in each point of space-time 
the time-like axis is the unit 4-velocity of the observer, while the spatial 
axes are a (gauge) convention for observer's gyroscopes. This framework was 
developed in the second and third paper of Refs.[17]. 

Even if the action of ADM tetrad gravity depends upon 16 fields, the 
counting of the physical degrees of freedom of the gravitational field does 
not change, because this action is invariant not only under the group of 
4-difcomorphisms but also under the 0(3,1) gauge group of the Newman- 
Penrose approach [68] (the extra gauge freedom acting on the tetrads in the 
tangent space of each point of space-time). 

The cotetrads E A a) (r,a r ) are the new configuration variables. They are 

o(«) 

connected to cotetrads A E A (t, a r ) adapted to the 3+1 splitting of space- 

oA 

time, namely such that the inverse adapted time-like tetrad E^(t, cr r ) is 
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the unit normal to the 3-space S T , by a standard Wigner boosts for time-like 
Poincare orbits with parameters <p^ (r, <r r ), a = 1, 2, 3 

a) (a) 4° (a) 4 4 ° W 

(v( a \z(cr)) + v c ' ) )(y w (z(cT)) + v w ) 

1 + V(°)(z(a)) ■ [ ' 

In each tangent plane to a point of S T this point-dependent standard Wigner 

o(«0 

boost sends the unit future-pointing time- like vector V = (1;0) into the 

unit time-like vector = 4 £^ a) l A = + f\ a y f (a) = -e ¥>(<*)) • As 

a consequence, the flat indices (a) of the adapted tetrads and cotetrads and 
of the triads and cotriads on S T transform as Wigner spin-1 indices under 
point-dependent SO(3) Wigner rotations Rr a \r b \(V(z(a)); A(z{a))) associ- 
ated with Lorentz transformations A^ a h^(z) in the tangent plane to the 
space-time in the given point of E T . Instead the index (o) of the adapted 
tetrads and cotetrads is a local Lorentz scalar index. 

The adapted tetrads and cotetrads have the expression 



4 



a 

= (l + n)(l;0) =cU, ^ = (ra( o) ; 3 e (a)r ), (23) 



where 3 e7 a ) and 3 e( a ) r are triads and cotriads on S T and n( ) = n r "e', a ^ 



3„r 

' % 

nT 3e (a)r 22 are adapted shift functions. In Eqs.(23) N(t, cr) = l+n(r, er) > 0, 
with n(r, <r) vanishing at spatial infinity (absence of super-translations), 
so that N(T,cr)dT is positive from S T to E T+ d T , is the lapse function; 
N r (r 7 cr) = n r (r, cr), vanishing at spatial infinity (absence of super-translations), 
are the shift functions. 

o o 

The adapted tetrads 4 -E( a ) are defined modulo SO(3) rotations 4 E^ — 

° A 

Ef, %)((,)(«( e ))%)> 3 e[ Q ) = Eb R (a)(b)(a(e)) 3 e r (by where a (a) (r,cr) are 



22 Since one uses the positive-definite 3-metric 5( Q )(M, one will use only lower flat spatial 
indices. Therefore for the cotriads one uses the notation 3 ei a ' d — 3 er„\ r with 5r, 



3e (a) \b)r- 



°{a)r °(a)(b) 
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three point-dependent Euler angles. After having chosen an arbitrary point- 
dependent origin ai a > (r, a) = 0, one arrives at the following adapted tetrads 
and cotetrads [n (a) = £ 6 n (b) i? (6)(o) (a (e) ) , Ea n (a) 3 e[ a) = Ea "(a) 3 e[ a) ] 



4 C) = 4 C = a; - E *(-) 3e 7a)) = ^ %«) = (°; 3e 7a)) 5 



+ n 

a 

o(o) (o) o(a) 

4 S A = i E A =(l + n)(l;0)=eI Al 4 £ A = (n (o) ; 3 e (a)r ), (24) 

which one will use as a reference standard. 
The expression for the general tetrad 



+ E 4 E {h) Rf b)ia) (a (c) )L( a \ a) ( nc) ), (25) 

afr 



shows that every point-dependent Lorentz transformation A in the tangent 
planes may be parametrized with the (Wigner) boost parameters <p( Q ) and 
the Euler angles a( a ), being the product A = RL of a rotation and a boost. 

The future-oriented unit normal to S T and the projector on S T are = 
e (1 + n) (l; o) , V s U Z B = e, ^ - e (1 + n) V T = ^ = 

T^(l;-E a S(a) S e[ a) ),M = ^-^Ai fl . 
The 4-metric has the following expression 

Vr = e [(1 + nf - V s n r n s ] = e [(1 + nf - £ n 2 o) ], 

a 

g rr = -en r = -e^2 n {a) 3 e(a)r, 

a 

3 x ^ 3 3 v 3_ 3_ 

9rs e 9rs e / J e (a)r e (a)s e / J e (a)r e (a)s) 

a a 

4 rr _ f 4 rr _ _ ^ _ _ Ea 3e (q) "(a) 

9 ~ (1 + n) 2 ' 5 6 (1 + n) 2 ~ 6 (1 + n) 2 ' 

4 ^ rs c(3r, rs ^ ^ \ * 3-r 3-s /r »(a) »(b) > 

9 - (1 + n)2 ) = - £ 2^ e (°) e (fc)( g WW ~ (i + n )2 )' 



4 



4 
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'N = -J=£= = V7(l + n) = 3 e(l + n), 
3 .9 = 7=( 3 e) 2 , 3 e = det 3 e (a)r . (26) 



The 3-metric 3 g rs has signature (+ + +), so that one may put all the flat 
3-indices down. One has 3 g ru3 g us = 5^. 



5.2 The ADM Phase Space and the ADM Hamilton 
Equations 

The given parametrization of the cotetrad fields leads to rewrite the ac- 
tion of ADM tetrad gravity in terms of the following 16 fields as configu- 
ration variables: three boost parameters <£( )(t, <j u ); the lapse N(t 7 <j u ) — 
1 + ti(t, a u ) and shift n( a )(r, <j u ) functions; the nine components of cotriad 
fields 3 e( a ) r (r, a u ) on the 3-spaces S T . As shown in the second and third pa- 
per of Ref.[17], in Rcf.[36] and in the first paper of Ref.[40], the ADM action 
for the gravitational field has the expression 



SgraV = 16VG / dTd3<7 V + ^ 3e£ ( a )( b )( c ) 3f2 rs(c) + 

+ 2(1 + n) ( 3g ° 1 )(°)(W°)('Q 3e w( n ('0l'- ~ dr 3e ( a )^ 
(d)( n (c)\s - d T 3 e( c ) s 



3 e s (d) (n (c)ls -d T 3 e (c)s ) (r,a u ). (27) 



In it 3 Q rs ( a ) = d r 3 uj s ( a) -d s 3 w r ( a ) -e( a )(6)(c) 3 ^r(b) 3 ^ s ( c ) is the field strength 
associated with the 3-spin connection 3 u! r r a \ = h^( a )(b)(c) 3e fb) 



'(6) \ Ur e (c)tt 



d u 3 e (c)r ) + 5 3 e^ fc) 3 e^ c) 3 e (d)r (d„ 3 e (d)u - d u 3 e (d)t ,) and ( 3 G D 1 )( a )(6)( c )(d) = 
^(a)(c)<5(6)(d) + <5(a)(d)<5(6)( c ) - 25( )((,)5( c ) (d ) is the flat (with lower indices) in- 
verse of the flat Wheeler-DeWitt super-metric 3 G ( a )(b)(c)(d) = ^(a)(c) 5(6)(d) + 

<*(a)(d) <*(6)(c) - 5(a)(6) *(c)(d), 3 G o(a ) (b)( e )(/) 3 G~ ( e)(/)( c )( d ) = 2 (<5 (a)(c) S (b){d) + 
5(a)(d) <*(6)(c))- 

The canonical momenta 71"^,^ (r, cr"), 7r„(r, ct"), 7r„ (a) (t, cr 11 ), 3 7t^(t, c"), 
conjugate to the configuration variables satisfy 14 first-class constraints: the 
ten primary constraints (the last three constraints generate rotations on quan- 
tities with fiat indices (a) like the cotriads) 



T <p w ( T >v U ) ~ °i n n (T,a u ) w 0, 7r„ (a) (r,CT u ) 
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3 M {a) (r,a u ) = e (a)mc) 3 e (b)r (T,a u ) 3 nl c) (T,a u ) w 0, (28) 
and the secondary supcr-Hamiltonian and super-momentum constraints 



n(r,a u ) 



n {a) (T,a u ) 



16ttG 
2ttG 



3„ , 3„r 3 s 3 o 

e £(o)(6)(c) e (a) e (6) J/ rs(c) 



3 G ( a )((,)( c )(d) 3 e(a)r ^(b) 3 e( c ) s 3 7T( d) 



(r,a")+X(r,a")«0, 



3 r 



3_r I 3 ^, r 



V ^(o) _ e (a)(6)(c) Wr(6) ^(c) + ^( a y 



e\ a) M r (r,O«0. 



(29) 



The functions .M(t, <t") and .A/f r (T, ct") describe the matter present in the 
space-time: A4(t, a u ) is the (matter- and metric-dependent) internal mass 
density, while A4 r (r, <r u ) is the universal (metric- independent) internal mo- 
mentum density. If the action of matter is added to Eq.(27), one can evaluate 



the energy-momentum tensor T ab (t,(j u ) 



2 8S„ 



(r, er) of the 



matter 23 and determine these functions from the following parametrization 

M(t, a u ) 



T tt (t, a u ) 
T Tr \r, a u ) 



3 e(l + n) 2 ]( T , a") 

(a) 



3 „r 
e (a) 



(1 + n) 3 ef } M s - n (a) M 



3 e(l + n) 2 



(T,a«). 



(30) 



The extrinsic curvature tensor of the 3-spaces S T as 3-manifolds embed- 
ded into the space-time has the following expression in terms of the barred 
cotriads of Eq.(24) and their conjugate barred momenta 



± V r c _ „ 



abu 



3 e(6) s + 3 e (a ) s 3 e (; ,) I ,^ 3 e (a)u : 



3 s 3 s 3- -u 
e(a)r e( a ) s e( fc )„7T (b) 



(31) 



Therefore the basis of canonical variables for this formulation of tetrad 
gravity, naturally adapted to 7 of the 14 first-class constraints, is 





n 


»(«) 


3 e (a)r 




TTn ~ 


7Tn (a) ~ 


3 7T r 



(32) 



23 The Hamilton equations imply 4 T AB = in accord with Einstein's equations and 
the Bianchi identity. 
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The behavior of these fields at spatial infinity (compatible with the absence 
of super-translations) is given in Eqs.(5.5) of the third paper in Refs.[17]; in 

particular for the cotriads one has 3 e( a ) r (r, a r ) -^ r ^oo ^1 + co1 ^ t ' ^ S ar + 
0(r" 3 / 2 ) (r = VLFF). 



From the action (29), after having added the matter action, one can obtain 
the standard non-Hamiltonian ADM equations (|r denotes the 3-covariant 
derivative in the 3-space S T with 3- metric 3 g rs ] 3 R rs is the 3-Ricci tensor of 



d T 3 g rs = n r \ s + n s \ r - 2 (1 + n) 3 K, 



d T 3 K rs = (1 + n) ( 3 R rs + 3 K 3 K, 



— n 



\s\r 



-.u 3 j 



K +n u i 3 K ^n ui K i 



2 K TU K s ) 



u 3 j 



(33) 



with the quantities appearing in these equations re-expressed in terms of the 
configurational variables of Eq.(32). 

Instead at the Hamiltonian level one can get the Hamilton equations for 
all the variables of the canonical basis (32), as shown in the first paper of 
Ref.[40], by using the Dirac Hamiltonian. As shown in Refs.[17], the Dirac 
Hamiltonian has the form (if the matter contains the electro-magnetic field 
there are extra terms with the electro-magnetic first-class constraints) 



1 , 

Hd = - Eadm 
c 



nH -n {a) n {a ) (t,ct") + 



+ 



K TTn + A„ (a) 7T„ (a) + \ V{a) TT V{a) + ^( a ) 3 M( a ) 



(t.O,(34) 



where Eadm is the weak ADM energy and the A's are arbitrary Dirac mul- 
tipliers. 



See Eqs. (2.22), (3.43)and (3.47) of the first paper of Ref . [40] for the expres- 
sion of the ten weak asymptotic ADM Poincare generators Eadm, Pad mi 
Jadmi ^"Adm- Since one is in a non-inertial rest frame (due to the absence 
of super-translations), one has the rest- frame conditions Padm ~ like in 
SR. Then one has to add the conditions K7 ADM « to eliminate the internal 
3-center of mass of the 3- universe like in SR [21]. Therefore the 3- universe 
can be seen as a decoupled external canonical non-covariant center of mass 
carrying a pole-dipole structure: the invariant mass Mc = i Eadm and the 
rest spin J a s dm ■ This view is in accord with an old suggestions of Dirac [8] . 
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In Ref.[40] the study of ADM canonical tetrad gravity was done with 
the following type of matter: N charged scalar particles (described by the 
canonical variables 7?[(t), K ir (r)) and the electro-magnetic field in the non- 
covariant radiation gauge (described by the canonical variables Aj_(r, cr u ), 
7r^_(r, a u ) as shown in Ref.[21]). The particles (described by an action like the 
one in Eq.(6)) have not only Grassmann- valued electric charges Qi (Qf = 0, 
Qi Qj — Qj Qi for i 7^ j) to regularize the electro-magnetic self-energies, but 
also Grassmann-valued signs of the energy (ijf — 0, m r/j = r]j m for i ^ j) to 
regularize the gravitational self-energies 24 . 

Instead in Ref.[41] the matter is a perfect fluid described by the action of 
Ref.[43] re-expressed in the 3+1 point of view in Rcfs.[42]. 

In the case of N particles the functions M. and M. r have the expression 
(see Ref. [40] for their form in presence of the electro-magnetic field) 



M(r, a") = £ 5\a\ (r)) m ^jm\ J + ^(r, a") « ir (r) ^ a) (r, a-) « ifl (r), 
»=i 

N 

M r (r, a u ) = ^Vi Kir(r), (35) 



i=i 



while in the case of dust [41], described by canonical coordinates a l (r, a u ), 
IIi(T } a u ), i = 1,2,3, they have the expression 



M(t, o u ) = [det (8 S ai)f + J2 Qa 2 Vra V sa d r a* d s ao n, n^r, a"), 

y arsij 

M r (T,a u )=J29 r a i (r,a u )n i (T,a u ). (36) 



Both quantities are two-valued. The elementary electric charges are Q = ±e, with e 
the electron charge. Analogously the sign of the energy of a particle is a topological two- 
valued number (the two branches of the mass-shell hypcrboloid) . The formal quantization 
of these Grassmann variables gives two-level fermionic oscillators. At the classical level 
the self-cncrgics make the classical equations of motion ill-defined on the world-lines of 
the particles. The ultraviolet and infrared Grassmann regularization allows to cure this 
problem and to get consistent solution of regularized equations of motion. Sec Refs.[28] for 
the electro-magnetic case. 
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6 The York Canonical Basis and the Inertial and Tidal 
Degrees of Freedom of the Gravitational Field 

The presence of 14 first-class constraints in the phase space having the 32 
fields of Eq.(32) as a canonical basis implies that there are 14 gauge variables 
describing inertial effects and 2 canonical pairs of physical degrees of freedom 
describing the tidal effects of the gravitational held (namely gravitational 
waves in the weak held limit). To disentangle the inertial effects from the 
tidal ones one needs a canonical transformation to a new canonical basis 
adapted to all the ten primary constraints (28) and containing the barred 
variables defined in Eq.(24). This is the topic of this Section. 



6.1 The York Canonical Basis 

A canonical transformation adapted to the ten primary constraints (28) was 
found in Ref.[36]. It implements the York map of Ref.[37] in the cases in 
which the 3-metric 3 g rs has three distinct eigenvalues and diagonalizes the 
York-Lichnerowicz approach (see Ref.[38] for a review). 

As said before Eq.(24), one can decompose the cotriads on S T in the 
product of a rotation matrix, belonging to the subgroup SO (3) of the tetrad 
gauge group and depending on three Euler angles (r, <r r ), and of barred 
cotriads depending only on six independent fields. The canonical transfor- 
mation Abelianizes the constraints 3 M( a )(r, a u ) w of Eqs.(28), satisfy- 
ing { 3 M {a) (T,a u ), 3 M {b) (T,a' u )} = C(o)(6)(c) 3 M (c) (r, a")<5 3 (a", a'"), and re- 
places them with the vanishing of the three momenta n^j | (r, a r ) w conju- 
gate to the Euler angles. 

The new canonical basis, named York canonical basis, is (a = 1, 2, 3; a = 
1.2) 



V(a) a (a) 


n n (a) 







Ra 


~ *($ ~ 


TTn ~ 7T fi(a) ss 






n-a 



(37) 

In it the cotriads and the components of the 4-metric have the following 
expression 

3e (a)r =J2 R (a)(b)(a {c )) 3 e { b)r - fi (a)(6) ("(c) ) Kb(^) ^ 7a ° R& , 

b b 
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V T = e[(l + n) 2 -^4 } ], 

a 

4 9rr = ~e ^ n (a) ^ e (a)r = ~ € ^2 ™( a ) ^^(^)r-> ^ = ^ = V det3 9rs, 

a a 

4 9rs = -e 3 g rs = -e4> 2/3 V ra {P)V sa {6 l )Q 2 a , Q a = e^' 2 ,(38) 

a 

The set of numerical parameters 7 5a appearing in Q a satisfies [17] E u -fan = 
0, E„ lau Tbu = hi, Ea 7s« lav = &uv ~ \ ■ Each solution of these equations 
defines a different York canonical basis. 

This canonical basis has been found due to the fact that the 3-metric 
3 g rs is a real symmetric 3x3 matrix, which may be diagonalized with an 
orthogonal matrix V(9 r ), V^ 1 = V T (E„ V ua Kfc = S ab , E a v ua V va = 5 UV , 
J2uv e wuv V ua V vb = E c e abc V cw ), detV — 1, depending on three parameters 
9 1 (i — 1,2,3) 25 , whose conjugate momenta Ilf^ are to be determined 
as solutions of the super-momentum constraints. If one chooses these three 
gauge parameters to be Euler angles #*(t, er), one gets a description of the 
3-coordinate systems on S T from a local point of view, because they give the 
orientation of the tangents to the three 3-coordinate lines through each point. 
However, for the calculations (see Refs.[40]) it is more convenient to choose 
the three gauge parameters as first kind coordinates 9 t (r 7 a) (—00 < 9 l < 
+00) on the 0(3) group manifold, so that by definition one has V ru (9 1 ) = 
( e - Hi T i 8 \ j where (Tj) ru — e rui are the generators of the o(3) Lie algebra 

V / ru 

in the adjoint representation, and the Euler angles may be expressed as 9 l = 
f l (9 n ). The Cartan matrix has the form A(9 n ) = S? ^ — and satisfies 
A ri {6 n )6 i = 5„9*; B{9 1 ) = A- 1 {9 1 ). 

From now on for the sake of notational simplicity the symbol V will mean 
V(9 l ). 

The extrinsic curvature tensor of the 3-space E T has the expression 



3 tf rs (r,<7") 



4ttG 



+ Y,Q"Qb (V ra V sb + V rb V sa ) 



(0) 



ab 



twi 



e abt Vwt Biw ?T 
Qb Qa — Qa Qfr 



(r,a u ). 
(39) 



25 Due to the positive signature of the 3-metric, one defines the matrix V with the fol- 
lowing indices: V ru . Since the choice of Shanmugadhasan canonical bases breaks manifest 
covariance, one will use the notation V ua = V uv <5„( a ) instead of V u ( a y 
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This canonical transformation realizes a York map because the gauge vari- 
able 7r^ (describing the freedom in the choice of the trace of the extrinsic 
curvature of the instantaneous 3-spaces S T ) is proportional to York internal 
extrinsic time 3 K. It is the only gauge variable among the momenta: this is a 
reflex of the Lorentz signature of space-time, because tt^ and n can be used 
as a set of 4-coordinates for the space-time [35] . The York time describes the 
effect of gauge transformations producing a deformation of the shape of the 
3-space along the 4-normal to the 3-space as a 3-sub-manifold of space-time. 

Its conjugate variable, to be determined by the super-Hamiltonian con- 
straint (interpreted as the Lichnerowicz equation), is <fr — (f> 6 — 3 e = 
^/det 3 g rs , which is proportional to Misner's internal intrinsic time; more- 
over <j) is the 3-volume density on S T : Vr = J r d 3 a <f>, R C S T . Since one has 
3 g rs = (j) 2 / 3 3 g rs with det 3 g rs — 1, (f> is also called the conformal factor of the 
3-metric. 

The two pairs of canonical variables i? a , n a , a — 1,2, describe the gen- 
eralized tidal effects, namely the independent physical degrees of freedom of 
the gravitational field. They are 3-scalars on S T and the configuration tidal 
variables i? a parametrize the two eigenvalues of the 3-metric 3 g rs with unit 
determinant. They are Dirac observables only with respect to the gauge trans- 
formations generated by 10 of the 14 first class constraints. Let us remark 
that, if one fixes completely the gauge and one goes to Dirac brackets, then 
the only surviving dynamical variables i? a and il a become two pairs of non 
canonical Dirac observables for that gauge: the two pairs of canonical Dirac 
observables have to be found as a Darboux basis of the copy of the reduced 
phase space identified by the gauge and they will be (in general non-local) 
functionals of the i? a , 7T a variables. 

Therefore, the 14 arbitrary gauge variables are ip^ (r, cr"), a( a )(r, cr 11 ), 
n(r, cr"), n( a )(r, cr 11 ), # 4 (r, cr"), tt^(t, cr 11 ): they describe the following general- 
ized inertial effects [36]: 

a) ct( a }{T, cr") and <p( a )(T, cr") are the 6 configuration variables parametriz- 
ing the 0(3,1) gauge freedom in the choice of the tetrads in the tangent plane 
to each point of S T and describe the arbitrariness in the choice of a tetrad 
to be associated to a time-like observer, whose world-line goes through the 
point (r, cr). They fix the unit ^-velocity of the observer and the conventions 
for the orientation of three gyroscopes and their transport along the world-line 
of the observer. The Schwinger time gauges are defined by the gauge fixings 
a {a )(r,o- u ) w 0, (p {a) (T,a u ) w 0. 

b) #*(t, cr") (depending only on the 3-metric) describe the arbitrariness in 
the choice of the 3-coordinates in the instantaneous 3-spaces E T of the chosen 
non-inertial frame centered on an arbitrary time-like observer. Their choice 
will induce a pattern of relativistic inertial forces for the gravitational field, 
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whose potentials are the functions V ra {9 1 ) present in the weak ADM energy 
Eadm- 

c) ni a \{r, <J U ), the shift functions, describe which points on different instan- 
taneous 3-spaces have the same numerical value of the 3-coordinates. They are 
the inertial potentials describing the effects of the non- vanishing off-diagonal 
components 4 <? rr (T, <j u ) of the 4-metric, namely they are the gravito-magnetic 
potentials 26 responsible of effects like the dragging of inertial frames (Lens- 
Thirring effect) in the post-Newtonian approximation. The shift functions 
are determined by the r-preservation of the gauge fixings determining the 
gauge variables 9 l (r, a u ) . 

d) TT^(T,a u ), i.e. the York time 3 K(t, <t u ), describes the non-dynamical 
arbitrariness in the choice of the convention for the synchronization of distant 
clocks which remains in the transition from SR to GR. Since the York time 
is present in the Dirac Hamiltonian, it is a new inertial potential connected 
to the problem of the relativistic freedom in the choice of the shape of the 
instantaneous 3-space, which has no Newtonian analogue (in Galilei space- 
time time is absolute and there is an absolute notion of Euclidean 3-space). 
Its effects are completely unexplored. Instead the other components of the 
extrinsic curvature of S T are dynamically determined once a 3-coordinate 
system has been chosen in the 3-space. 

e) I + n(r, a u ), the lapse function appearing in the Dirac Hamiltonian, 
describes the arbitrariness in the choice of the unit of proper time in each 
point of the simultaneity surfaces S T , namely how these surfaces are packed 
in the 3+1 splitting. The lapse function is determined by the T-preservation 
of the gauge fixing for the gauge variable 3 K (r, a u ). 

As shown in Ref.[35], the dynamical nature of space-time implies that 
each solution (i.e. an Einstein 4-geometry) of Einstein's equations (or of the 
associated ADM Hamilton equations) dynamically selects a preferred 3+1 
splitting of the space-time, namely in GR the instantaneous 3-spaces are 
dynamically determined modulo only one inertial gauge function (the gauge 
freedom in clock synchronization in GR). In the York canonical basis this 
function is the York time, namely the trace of the extrinsic curvature of the 
3-space. Instead in SR the gauge freedom in clock synchronization depends on 
four basic gauge functions, the embeddings z M (r, er r ), and both the 4-mctric 
and the whole extrinsic curvature tensor were derived inertial potentials. 
Instead in GR the extrinsic curvature tensor of the 3-spaces is a mixture 



6 In the post-Newtonian approximation in harmonic gauges they are the counterpart 
of the electro-magnetic vector potentials describing magnetic fields [38]: A) N = 1 + n, 

n d = — with <&q the gravito- electric potential; B) n r d = ^ Aq t with Aq t the 

gravito-magnetic potential; C) Eq t = d r <Pg ~ 9 T (5 Aq t ) (the gravito- electric field) and 
Bq r = e ruv d u Aq v = c Qq r (the gravito-magnetic field). Let us remark that in arbitrary 
gauges the analogy with electro-magnetism breaks down. 
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of dynamical (tidal) pieces and inertial gauge variables playing the role of 
inertial potentials. 



6.2 3- Orthogonal Schwinger Time Gauges and 
Hamilton Equations 

As shown in the first paper in Rcfs.[40], in the York canonical basis the Dirac 
Hamiltonian (34) becomes (the A's are arbitrary Dirac multipliers; the Dirac 
multiplier A r (r) implements the rest frame condition P T ADM ~ 0) 



H D = - Eadm + I d 4 a 



nU - n (a) U {a) (r, a w ) + \ r {r) P r ADM + 



d A a 



(r,0, (40) 



with the following expression for the weak ADM energy 



E 



ADM 



6 



16ttG 



+ 2 E 

abtwiuvj 



^abt tabu Vwt ^iw Vvu Bjv 



) ](r.O- 



(41) 



In it 5(r, a") is a function of 0, 0* and i? a (given in Eq.(B8) of the first 
paper in Ref. [40]), which play the role of an inertial potential depending on 
the choice of the 3-coordinates in the 3-space (it is the r — r term in the 
scalar 3-curvature of the 3-space). 

Eq.(41) shows that the kinetic term, quadratic in the momenta, is not 
positive definite. While the kinetic energy of the tidal variables and the last 
term 27 are positive definite, there is the negative kinetic terms (vanishing 



only in the gauges 3 K(t, a u ) =0) 



24ttG 



/ D 3 a<j>(T,o- u ) 3 K 2 (T,o- u ). It is an 



inertial potential associated with the inertial gauge variable York time, which 
is a momentum due to the Lorentz signature of space-time. It was known that 
this quadratic form is not definite positive, but only in the York canonical 
basis this can be made explicit. 



In the York canonical basis it is possible to follow the procedure for the 
fixation of a gauge natural from the point of view of constraint theory when 
there are chains of first-class constraints [9]. This procedure implies that 



27 It describes gravito-magnetic effects. 
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one has to add six gauge fixings to the primary constraints without without 
secondaries (itip^ (t, <j u ) rs 0, 7r Q(a) (t, cr u ) w 0) and four gauge fixings to 
the secondary super-Hamiltonian and super-momentum constraints. These 
ten gauge fixings must be preserved in time, namely their Poisson brackets 
with the Dirac Hamiltonian must vanish. The T-preservation of the six gauge 
fixings determining the gauge variables a^(r, a u ) and <P( )(r, cr u ) produces 
the equations determining the six Dirac multipliers A v (r, a u ), A a , o) (r, a u ). 
The r-preservation of the other four gauge fixings, determining the gauge 
variables 9 1 (t,<t u ) and the York time 3 K(r,a u ), produces four secondary 
gauge fixing constraints for the determination of the lapse and shift functions. 
Then the T-preservation of these secondary gauge fixings determines the four 
Dirac multipliers A„(r, a u ), A S(a) (r, a u ). Instead in numerical gravity one 
gives independent gauge fixings for both the primary and secondary gauge 
variables in such a way to minimize the computer time. 

In Section V of the first paper in Refs.[40] there is a review of the gauges 
usually used in canonical gravity. It is shown that the commonly used family 
of the harmonic gauges is not natural according to the above procedure. 
The harmonic gauge fixings imply hyperbolic PDE for the lapse and shift 
functions, to be added to the hyperbolic PDE for the tidal variables. Therefore 
in harmonic gauges both the tidal variables and the lapse and shift functions 
depend (in a retarded way) from the no-incoming radiation condition on the 
Cauchy surface in the past (so that the knowledge of 3 K from the initial time 
till today is needed). 

Instead the natural gauge fixings in the York canonical basis of ADM 
tetrad gravity are the family of Schwinger time gauges, where the 0(3,1) 
gauge freedom of the tetrads is eliminated with the gauge fixings (implying 

A ¥ , (o) (T,<7«) = A a(o) (r,a«) = 0) 

«(.)(r,ff*)«0, ^ (a) (r,O«0, (42) 
and the subfamily of the 3-orthogonal gauges 

6\t, a u ) w 0, 3 K(t, a 11 ) w F(t, a u ) = numerical function, (43) 

in which the 3-coordinates are chosen in such a way the the 3-metric in the 
3-spaccs E T is diagonal. The r-preservation of Eqs.(43) gives four coupled 
elliptic PDE for the lapse and shift functions. Therefore in these gauges only 
the tidal variables (the gravitational waves after linearization) , and therefore 
only the eigenvalues of the 3-metric with unit determinant inside E T> depend 
(in a retarded way) on the no-incoming radiation condition. The solutions <j> 
and 7r 4 of the constraints and the lapse 1+n and shift n< a \ functions depend 
only on the 3-space U T with fixed r. If the matter consists of positive energy 
particles (with a Grassmann regularization of the gravitational self-energies) 
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[40] these solutions will contain action-at-a-distance gravitational potentials 
(replacing the Newton ones) and gravito-magnetic potentials. 

In the family of 3-orthogonal gauges the weak ADM energy and the super- 
Hamiltonian and super-momentum constraints (they are coupled elliptic PDE 
for their unknowns) have the expression (see Eq.(3.47) of the first paper in 
Ref . [40] for the other weak ADM Poincare generators) 



EadmW=q = c J d a M\oi= 
2irG - ! 



16ttG 



S|fl«=0 + 



+ ^<r 1 (-3(H) 2 + 2 E^ + 

b 

(6) (8) 

+ 2E _e^^ )](r/)i 



abij ^Q a Q b 1 — QbQa 1 



H(r,a u )\ ei=0 



+ 



16ttG 
2ttG - ! 



1/6 (t, a") [8 A ~4> 1,& ~ 3 ^ =0 4> 1/6 }(r, a") + M\ s , =0 (r, a u ) + 



+ 2 E 



abi i \Qa 1 — QbQa 1 



3N 2 + 2^U S 2 + 



^abi ^abj TTj 



W( o) |ji = o(r,a«) = r 2 (r,<r)[££ 



+ 2 EE 



QbQa 1 — QaQ 
£abi Q a 



b^a i {QbQa 1 - QaQb 1 ) 

+ Qa 1 (/ da^ + Y, (75a ^ 27 s - 9 a J? 5 n h ) + M a )] (r, a"), 



A = E 2 + 2 E 7ar 9 r J?a(T, (T u ) d r ' , 
r a 

s*=o(r, a-) = (>/ 3 £ 2 [jj 1 9 a 0) 2 + 

a 

+ E ( E ( 2 Tba Tea - he) 9a R~c ~ \ Tba 4>'' 9 a 4>) d a R b ] ) (t, <7 U ) . 
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(44) 

In the first paper in Refs.[40] there is the explicit form of the Hamilton 
equations for all the canonical variables of the gravitational field and of the 
matter replacing the standard 12 ADM equations and the matter equations 
4 VaT ab = in the Schwingcr time gauges and their restriction to the 
3-orthogonal gauges. They could also be obtained from the effective Dirac 
Hamiltonian of the 3-orthogonal gauges, which is evaluated by means of a r- 
dependent canonical transformation sending the gauge momentum 7t^(t, <j u ) 

in the gauge-fixing conditions tt~(t, a u ) = ^ G (^K(r,a u ) — F(t, a u )^ « 
and which is given in Eq.(4.39) of the second paper of Ref.[40]. 

These equations are divided in five groups: 

A) The contracted Bianchi identities, namely the evolution equations for 
the solutions 0(t, cr") and tt\ (t,<t u ) of the super-Hamiltonian and super- 
momentum constraints: they are identities saying that, given a solution of 
the constraints on a Cauchy surface, it remains a solution also at later times. 

B) The evolution equation for the four basic gauge variables 1 (t, a u ) and 
3 K(t, a u ) (the equation for the York time is the Raychaudhuri equation 28 ): 
these equations determine the lapse and the shift functions once four gauge- 
fixings for the basic gauge variables are given. 

C) The equations d T n(r, a u ) = A„(r, <r u ) and d T n( a ){ T i °~ u ) = ^n (a) (t, <t u ). 
Once the lapse and shift functions of the chosen gauge have been found, they 
determine the associated Dirac multipliers. 

D) The hyperbolic evolution PDE for the tidal variables i? a (T, a u ), n a (T, a u ). 
When the equations for d T R a (T 7 a u ) is inverted to get i7 a (T, a u ) in terms 
of Ra{T,a u ) and its derivatives, then the Hamilton equations for -Z7 5 (t, a u ) 
become hyperbolic PDE for the evolution of the physical tidal variable 
Ra(r,a u ). 

E) The Hamilton equations for matter, when present. 

Given a solution of the super-momentum and super-Hamiltonian con- 
straints and the Cauchy data for the tidal variables on an initial 3-space, one 
can find a solution of Einstein's equations in radar 4-coordinates adapted to 
a time-like observer in the chosen gauge. 



28 This equation is relevant for studying the developments of caustics in a congruence of 
time-like geodesies for converging values of the expansion 8 and of singularities in Einstein 
space-times [69]. However the boundary conditions of asymptotically Minkowskian space- 
times without super-translations should avoid the singularity theorems as it happens with 
their subfamily without matter of Ref.[16]. 
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6.3 The Congruence of Eulerian Observers and the 
non-Hamiltonian First- Order ADM Equations of 
Cosmological Spacetimes 

Like in SR one can consider the congruence of the Eulerian observers with zero 
vorticity associated with the 3+1 splitting of space-time, whose properties 
are described by Eq.(5). In the first paper of Ref.[40] it is shown that in 
ADM tetrad gravity the congruence has the following properties in each point 
(r;^) 29 : 

a) The acceleration 3 a A = 1 B4 VbI A = 4 g AB3 aB has the components 
3 a T = 0, V = e<T 2 / 3 Q- 2 V ra V sa d s ln{\ + n), V = Q~ x V ra n {a) 
d r ln(l + n), 3 a r = —d r In (1 + n). 

b) The expansion 30 coincides with the York time: 

e = 4 V A l A = -e*K = -e 1 -^^. (45) 

In cosmology the expansion is proportional to the Hubble constant and the 
dimcnsionless cosmological deceleration parameter is q = 3 l A 4 V+ g — 1 = 
-39- 2 l A d A 9-l. 

c) By using Eqs.(24) it can be shown that the shear 31 gab — &ba — 

o(«) o(/3) 

-f ( 3 a A Ib + 3 a B I A ) + § ( 4 V A h + 4 V B l A ) - \6 3 h AB = a {a)m 4 E A 4 E B 
has the following components cr (o)(o ) = cr (o)(r) = 0, cr( a )(6) = °~(b)(a) = 
( 3 K rs ~l 3 g rs 3 K) 3 e\ a) 3 e s {b) , ^ a CT (o)(o) = 0. o-( a )(b) depends upon the canon- 
ical variables 9 r , (f>, R a , and 7T a . 

By using Eqs.(39) for the extrinsic curvature tensor one finds that the 
diagonal elements C(a)(a) of the shear are also connected with the tidal mo- 
menta 77 a , while the non-diagonal elements C( a )(b) \a^b are connected with 
the momenta ir\ (the unknowns in the super-momentum constraints) 



c 3 



29 See the " 1+3 point of view" of Rcf. [70] for a discussion of gravity in terms of the second 
non-surface- forming congruence of time-like observers associated with a 3+1 splitting of 
space-time. 

30 It measures the average expansion of the infinitesimally nearby world-lines surrounding 
a given world-line in the congruence. 

31 It measures how an initial sphere in the tangent space to the given world-line, which is 
Lie-transported along the world-line tangent Z M (i.e. it has zero Lie derivative with respect 
to Z M c9 M ) , is distorted towards an ellipsoid with principal axes given by the eigenvectors of 
(T M „, with rate given by the eigenvalues of <r M „. 
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] K rs = 4> 2/3 J2(-^6 5 ab + cr {a){b ^Q a Q b v ra v sb ^ 

ab 

+ "(a)(6))- (46) 



>0 2 ^ 3 Qr V - 



Therefore the Eulerian observers associated to the 3+1 splitting of space- 
time induce a geometrical interpretation of some of the momenta of the York 
canonical basis: 

1) their expansion 8 is the gauge variable York time 3 K = 12 Js, de- 
termining the non- dynamical gauge part of the shape of the instantaneous 
3-spaces S T as a sub-manifold of space-time; 

2) the diagonal elements of their shear describe the tidal momenta Tla, 
while the non-diagonal elements are connected to the variables Trf\ deter- 
mined by the super-momentum constraints. 

In Eq.(44), valid in the 3-orthogonal gauges, the term quadratic in the 
momenta irf^ in the weak ADM energy and in the super-Hamiltonian con- 
straint can be written as j§^q 4* J2 a b a^b °fa)(6)> wnuc the super- momentum 
constraints can be written in the form of PDE for the non-diagonal elements 
of the shear 



b=£a 



db<J(a)(b) 



+ 1 + (Tba ~ lib) 9b Rbj V(a)(b) 

b 

~ ^ [* da^ + Y, (Tba da U g - d a R~ b 77,) + 

b 

+ M a ]){r,a u ) wO. (47) 



As a consequence, by using 3 g rs of Eq.(38) and 3 K rs of Eq.(46), the first- 
order non-Hamiltonian ADM equations (33) can be re-expressed in terms of 
the configurational variables n, ni a \, <j>, l , 7? a , and of the expansion 9 and 
shear o-< a \(h\ of the Eulerian observers. Then the 12 equations can be put in 
the form of equations determining d T <j>, d T 7? s , d T l , drO and 9 T o'(a)(6)- I n 
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Eqs.(2.17) of the first paper in Ref.[40] this manipulation is explicitly done 
for the first six equations (33). 

These results are important for extending the identification of the inertial 
and tidal variables of the gravitational field, achieved with the York canonical 
basis, to cosmological space-times. Since these space-times are only confor- 
mally asymptotically flat, the Hamiltonian formalism is not defined. However, 
they are globally hyperbolic and admit 3+1 splittings with the associated con- 
gruence of Eulerian observers. As a consequence, in them Einstein's equations 
arc usually replaced with the non-Hamiltonian first-order ADM equations 
plus the super-Hamiltonian and super-momentum constraints. Our analysis 
implies that, since the 4- metric can always be put in the form of Eqs.(38), 
the inertial gauge variables of the cosmological space-times are n, ni a -\ , 9 l and 
the expansion 6 = —e 3 K, while the physical tidal variables are R a and the 
diagonal components of the shear <J( a )( a ) (J2a <J (a.)(a) — 0). The unknown in 
the super-Hamiltonian constraint is the conformal factor <f> of the 3-metric in 
£ T , while the unknowns in the super-Hamiltonian constraints are the non- 
diagonal components of the shear C( )(6) U^6- 
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7 Post-Minkowskian Linearization in Non-Harmonic 
3-Orthogonal Gauges and Post-Minkowskian 
Gravitational Waves 

In the second paper of Ref.[40] it was shown that in the family of non- 
harmonic 3-orthogonal Schwinger gauges it is possible to define a consistent 
linearization of ADM canonical tetrad gravity plus matter (N charged scalar 
particles of masses m i} Grassmann- valued signs of energy Grassmann- 
valued electric charges Qi, plus the electro-magnetic held in the radiation 
gauge) in the weak held approximation, to obtain a formulation of Hamil- 
tonian Post-Minkowskian (HPM) gravity with non-flat Riemannian 3- spaces 
and asymptotic Minkowski background. 

In the standard linearization one introduces a fixed Minkowski background 
space-time, introduces the decomposition 4 g^ v (x) — 4 rj fl „ + 4 h^(x) in an 
incrtial frame and studies the linearized equations of motion for the small 
Minkowskian fields 4 /i^„(x). The approximation is assumed valid over a big 
enough characteristic length L interpretable as the reduced wavelength X/2tt 
of the resulting gravitational waves (GW) (only for distances higher of L the 
linearization breaks down and curved space-time effects become relevant). For 
the Solar System there is a PN approximation in harmonic gauges, which is 
adopted in the BCRS [5] and whose 3-spaces ts = const, have deviations of 
order c~ 2 from Euclidean 3-spaces. 

See Refs.[44, 71] and Appendix A of the second paper in Refs.[40] for a 
review of all the results of the standard approach and of the existing points 
of view on the subject [72, 73, 74, 75, 76]. 

In the class of asymptotically Minkowskian space-times without super- 
translations the 4-metric tends to an asymptotic Minkowski metric at spatial 
inhnity, gAB — > Vab, which can be used as an asymptotic background. 
The decomposition 4 gAB = 4 i]ab + 4 ^(i)ab, with a hrst order perturbation 
4 ^(i)AB vanishing at spatial infinity, is defined in a global non-inertial rest 
frame of an asymptotically Minkowskian space-time deviating for first order 
effects from a global inertial rest frame of an abstract Minkowski space-time 
-^(oo) • The non-Euclidean 3-spaces S T will deviate by first order effects from 
the Euclidean 3-spaces £ T (oo) of the inertial rest frame of coinciding 
with the limit of S T at spatial inhnity. When needed differential operators 
like the Laplacian in S T will be approximated with the flat Laplacian in 

If £ << 1 is a small a-dimensional parameter, a consistent Hamiltonian 
linearization implies the following restrictions on the variables of the York 
canonical basis in the family of 3-orthogonal gauges with 3 K(t, cr) = F(t, cr) 
= numerical function (in this Section one uses the notation cr for the curvi- 
linear 3-coordinates a r ) 
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R- a {T,cr) = R {1) - a {T,cr) = 0{() « 1, 
n- a {T,<r)=n {1) - a {T,a) = ^0{Q, 



4>(t, cr) = ^det 3 g rs (r, cr) = 1 + 6 (r, cr) + 0(C 2 ), 
JV(t, cr) = 1 + n(T, cr) = 1 + n ( i)(T, cr) + 0(C 2 ), 
e 4 9T r(r,cr) = l + e 4 ^ (1)Tr (r,cr) = 1 + 2 n (1) (r, cr) + 0(C 2 ), 
n (a) (r,cr) = -e i g Ta (T,(r) = -e 4 /i (1)Tr (r, cr) = n (1)(a) (r, cr) + 0(( 2 ), 



3 if (T, cr) = — — 7T^ ( T) cr) = 3 K (1) ( T) cr) = — 



(l)^(T,Cr) = z O(C) 



o-(o)(6)U#6(T,cr) = a- (1 )( a )( 6) | a ^ 6 (T,Cr) = y 0(C), 




e 5rs(T, cr) = c> rs - e /i(i) rs (r, <r) = 



= [1 + 2 (r« (r, cr) + 2 (r, cr))] <5 rs + 0(( 2 ), 

2 3 
ri 1 )(r,Cr)= TaaRair,*), R*(t, Ct) = £ 7 a„ T« (t, <t). (48) 



The tidal variables Ra(T,cr) are slowly varying over the length i and 
times L/c; one has (t^) 2 = 0(C)) where 4 7£ is the mean radius of curvature 
of space-time. 

The consistency of the Hamiltonian linearization requires the introduction 
of a ultra-violet cutoff M for matter. For the particles, described by the canon- 
ical variables ^(t) and ACj(r), this implies the conditions = 0(C)- 
With similar restrictions on the electro-magnetic field one gets that the 
energy-momentum tensor of matter is T ab (t, cr) = T^ b (t, er)+0(C 2 ). There- 
fore also the mass and momentum densities have the behavior M(r,cr) = 
M(i)(t, cr) + 0(C 2 ), M r (r, cr) = Al(i) r (r, cr) + 0(C 2 )- This approximation is 
not reliable at distances from the point particles less than the gravitational 
radius R M = ^ « 10~ 29 M determined by the cutoff mass. The weak ADM 
Poincare generators become equal to the Poincare generators of this matter 
in the inertia! rest frame of the Minkowski space-time plus terms of 
order 0(C 2 ) containing GW and matter. Finally the GW described by this 
linearization must have wavelengths satisfying X/2tt ps L >> Rm- If all the 
particles are contained in a compact set of radius l c (the source) , one must 
have l c » Rm for particles with relativistic velocities and l c > Rm for slow 
particles (like in binaries). See Ref.[44] for more details. 



ar— 1 
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With this Hamiltonian linearization one can avoid to make PN expansions: 
one gets fully relativistic expressions, i.e. a HPM formulation of gravity. 

The effective Hamiltonian adapted to the 3-orthogonal gauges and replac- 
ing the weak ADM energy is i (^Eadm(i)+Eadm(2)^J+j§^q J d 3 e (d T 3 K {1) ' 

1 (4 E Q d l r a 1] ~ *r ^(i))]) ( T ' ° u ) + 0(( 3 ) in the PM linearized theory. 

In the second paper of Refs.[40] one has found the solutions of the super- 
momentum and super-Hamiltonian constraints and of the equations for the 
lapse and shift functions with the Bianchi identities satisfied. Therefore 
one knows the first order quantities tt)iL(t, cr), </>(t, cr) = 1 + 6 4>(i)(t, cr), 
l + n(i)(r, cr), n(i)( a )(r, cr) (the quantities containing the action- at- a- distance 
part of the gravitational interaction in the 3-orthogonal gauges) with an ex- 
plicit expression for the PM Newton and gravito-magnetic potentials. In ab- 
sence of the electro-magnetic field they arc (the terms in (r, cr) describe 
the contribution of GW) 



1+ 



32 



4>{t,<t) = 1 + 6 ( i)(T,er) = 



1 + ^mf c 2 + ^(t) 

c 3 ^ * |<T-r7 4 (r)| 

3 /-^E^L^^en) 



8?r J |<7 — <Ti| 

e 4 3rr (t, cr) = 1 + 2 n (1) (t, cr) = 1 - 2 <9 T 3 /C (1) (r, cr) 



2G 



4-^ k-^(r)| l i+ mf C 2 + K |J' 



-e 4 5Ta( T , cO = n (1)(a) (t, cr) = a a 3 /C ( i)(r, cr) - 

lK-^(T))^(r).(cr-r, 4 (r))' 



/ 



d 3 a 1 



4ir cr — cr! 



3la <9 T 



2ri 1 )(r,«r 1 )- 



32 Quantities like \ri i (r) — rfj(r)\ are the Euclidean 3-distance between the two particles 
in the asymptotic 3-space S T ^ 00 ^ ) , which differs by quantities of order 0(Q from the real 
non-Euclidean 3-distance in E T as shown in Eq.(3.3) of the third paper in Ref.[40]. 
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/ 



« Ec^c^c (1) (r,«T 2 ) 
a a 2 



8ir \(Ti - <J2\ 

V(l)(a)(b)\a=£b( T > (T ) = 2 ( 9 <» "(!)(*>) + 5 fc "(1) (a) ) I a^b (t, Cr) . (49) 

Instead the linearization of the Hamilton equations for the tidal variables 
Ra(r, er) implies that they satisfy the following wave equation 33 (A and □ 
are the flat Laplacian and the flat D'Alambertian on S T ^ 00 ^) 



dlRa{T,(r) = AR a (T,er)+^2-faa d T d a n {l){a 



8ttG 



+ dl n (1) + 2 flg - 2 a Q 2 r« + — 2^ (r, <r). (50) 



By using Eqs.(49) this wave equation becomes 



M- a - b R- b (r,(T) = E a (T,<r), 

b 

^ "fba — 2 E 



M ab = S ab 



lib 



9 1 

A)' 



*m = ^ E ^ ( 4A4 w - 1 E »^d)c) + 



"cL d, 



9. 



□ ^M ab r b (1) (r,,T) = Y^TaaE- a (T,<r), 

b a 

lab = TaaTbbM a - b = S ab (l - 2 ^) + ^ (l , ^ , ^ 

^ M ab = 0, M s6 = ^ 7aa Tbb M ab . (51) 



u a \ °b 



ab 



33 For the tidal momenta one gets 77a(r, it) = [9 T -Ra - X] a ">' aa ^ a n (l)(a)]( T > <T ) + 

0(£ 2 ), so that the diagonal elements of the shear arc C(i)( a )( a ) (t, it) = [— 7a a d T Ra + 
"(l)(o) - 3 Eb ™(i)(i>)K T .' T ) +°(C 2 )- 
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To understand the meaning of the spatial operators M ab and M ab , one 



must consider the perturbation 4 /i(i) rs (T, cr) = —2e8 rs (r, 



(i) 



i))(t,ct) 



of Eq. (48) and apply to it the following decomposition, given in Ref [45] , 

4 fyl)rs(T , Cr) = ( 4 /lflJrs + g $rs #(1) + ^ i®r e (l)s + 9 S C(i) r ) + 

+ (9 r d s - ij„A)A ( i))(T,<T), (52) 

with ^ r <9 r £(!) r = and 4 hJ^ rs traceless and transverse (TT), i.e. ^ r 4 hJ^ rr 
O^Er d r 4 hf^ rs = 0. Since one finds (r, a ) = -12 e (1) (t, cr), A ( i)(t,ct) = 

"3e E„ fr^W) and e (1)r (r,cr) = -4e^ (r« - £„ § ^(r.ff), 
it turns out that the TT part of the spatial metric is independent from <f>^ 
and has the expression 



2^(rW + rV) + ^J2%rW}(T,*), 



A 



4 h TT 



(i)aa(r,cr) = -2eJ2M ab rl 1 \T,cr). 



(53) 



Therefore the spatial operator M ab connects the tidal variables -R a ( T , cr) of 
the York canonical basis to the TT components of the 3-metric. By applying 
the decomposition (52) to the spatial part T^(t, cr) of the energy-momentum 
one verifies that like in the harmonic gauges [44] the TT part of the 3-metric 



satisfies the wave equation □ 4 /i^ s T (r, cr) 



16ttG rp(TT) 
-(l)rs 



(r,cr). 



The retarded solution of the wave equation with a no-incoming radiation 
condition gives the following expression for the tidal variables (the HPM-GW) 



R- a {T, cr) = -J2 Taa ^ (r, cr)^ Taa (r, *) 

a 

2G 



8ttG 



/ 



T,\, (t- \<t - CT^CTi) 



(1) 



<T — <Ti 



n- a (r,cr) = (Y,M a - b d T R h -Y, 



4irG 1 



c 3 A 



t(49«M ( i)„ 



-^^a c X (1)c ) + a o 23 /C (1) ])(r,eT). 



(54) 
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The explicit form of the inverse operator is given in the second paper of 
Ref.[40]. By using the multipolar expansion of the energy- momentum 
of Ref.[58] in the HPM version adapted to the rest- frame instant form of 
dynamics of Ref. [26] , one gets 



G - d 2 g( TT ) aa l Tr (r — \tj\) 
Ra(T,cr) = ~— ^ JaaM^ 1 — : — : h (higher multipoles) , 



(55) 



where q( TT ) aa l TT (r) is the TT mass quadrupole with respect to the center 
of energy (put in the origin of the radar 4-coordinates). An analogous re- 
sult holds for 4 /i^ s t (t, <t) and this implies a HPM relativistic version of the 
standard mass quadrupole emission formula. 

Moreover, notwithstanding there is no gravitational self-energy due to the 
Grassmann regularization, the energy, 3-momentum and angular momentum 
balance equations in HPM-GW emission are verified by using the conserva- 
tion of the asymptotic ADM Poincare generators (the same happens with 
the asymptotic Larmor formula of the electro-magnetic case with Grassmann 
regularization as shown in the last paper of Ref. [27]). See Refs.[44, 76, 77] for 
the use of the self-energy in the standard derivation of this result by means 
of PN expansions. 

Eqs.(49) and (54) show that the HPM linearization with no- incoming radi- 
ation condition and Grassmann regularization is a theory with only dynamical 
matter interacting through suitable action-at-a-distance and retarded effec- 
tive potentials. Instead in relativistic atomic physics in SR the no-incoming 
radiation condition and the Grassmann regularization kill also the retarda- 
tion leaving only the action-at-a-distance inter-particle Coulomb plus Darwin 
potentials. See Eq.(7.22) of the second paper of Ref. [40] for the expression of 
the weak ADM energy till order 0(C 3 )- 

Moreover it can be shown that the coordinate transformation f = r, a r = 
°' r + \ TsT — J2 C % ^c 1 ^ (t, <t), introducing new r-dependent radar 

3-coordinates on the 3-space ZJ T , allows one to make a transition from the 3- 
orthogonal gauge with the 4-metric given by Eqs.(48) and (49) to a generalized 
non- 3 -orthogonal TT gauge containing the TT 3-metric (53) 



A 9(i)AB = 4 VAB + 

( -2 %- 3 ^ (1) + a(matter) -^ 3 K {1) + A r (r^ ] ) + ^(matter) 

+e 

\-f 3 ^ ( i) + A s (r^) + (3 s {matter) [^(rf 1 )) + 7 (matter)] 5 rs 



+ 
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+ o(c 2 ), 



-2 %- 3 K (1) + a(matter) 3 K (1) + ^(matter) \ 



4 3ab = 4 1iB + e | | + 

% 3 if ( i) + ft (matter) e 4 hf^ r s + S rs j (matter) ^ 

+ 0(C 2 ). (56) 



The functions appearing in Eqs.(56) are: A r (r^) = —\d T |f ^ 
£ (Ad + Ec T ( T))' ^(matter) = £ (4M (1)r - ff E c »^(i)o) 



8ttG 



A 

-{(matter) = ^G^M (1) . 



Also in absence of matter this TT gauge differs from the usual harmonic 
ones for the non-spatial terms depending upon the inertial gauge variable 
non-local York time 

3 IC w (r,a) = ^ s K {1) (r,a), (57) 

describing the HPM form of the gauge freedom in clock synchronization. 

If one uses the coordinate system of the generalized TT gauge, one can 
introduce the standard polarization pattern of GW for 4 /^T s T (see Refs.[44, 
45, 78]) and then the inverse transformation gives the polarization pattern of 
HPM-GW in the family of 3-orthogonal gauges. 

If the matter sources have a compact support and if the matter terms 
^ A / J(i)(t, <t) and -Mm r (T, a) are negligible in the radiation zone far away 
from the sources, then Eq.(56) gives a spatial TT-gauge with still the explicit 
dependence on the inertial gauge variable 3 K.^(t, er) (non existing in New- 
tonian gravity), which determines the non-Euclidean nature of the instanta- 
neous 3-spaces. Then one can study the far field of compact matter sources: 
the restriction to the Solar System of the resulting HPM 4-metric 34 is com- 
patible with the harmonic PN 4-metric of BCRS [5] if the non-local York 
time is of order c -2 . The resulting shift function should be used for the HPM 
description of gravito-magnetism (see Refs.[38, 79, 80] for the Lense-Thirring 
and other associated effects). 



34 See Eq.(7.20) of the second paper in Ref.[40], where 4 <? tt (t,<t) and 4 gTr(T,cr) are 
explicitly depending on the non-local York time. 
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The TT gauge allows one to reproduce the various descriptions of the 
GW detectors and of the reference frames used in GW detection in terms of 
HPM-GW: this is done in Subsection VIID of the second paper of Ref.[40], 
where the effect of a HPM-GW on a test mass is given in terms of the proper 
distance between two nearby geodesies. 

The HPM-GW propagate in non-Euclidean instantaneous 3-spaces S T dif- 
fering from the inertial asymptotic Euclidean 3-spaces E T (oo) a t the first 
order. In the family of 3-orthogonal gauges with York time 3 ^(i) (t, er) w 
Fm (t, er) = numerical function, the dynamically determined 3-spaces S T 

have an intrinsic 3-curvature 3 R\ei = o — 2 d% determined only by the 
HPM-GW (and therefore by the matter energy-momentum tensor in the past 
as shown by Eq.(54)). Their extrinsic curvature tensor as sub- manifolds of 
space-time is 

3 K(i) rs « cr(i)( r )( s )| r ^ s + S rs - d T -T r (1) + d r ri(i)( r ) - d a ri(i)(a)), 

a 

(58) 

with ri(i)( r ), cr(i)(r)(s)|r^s an d rj: given in Eqs. (49) and (54). The York time 
appears only in Eq.(58): all the other PM quantities depend on the non-local 
York time 3 JC (1) (r, er) w i (r, er) 

In the third paper of Refs.[40], where the matter is restricted only to the 
particles 35 , one evaluates all the properties of these HPM space-times: 

a) the 3- volume element, the 3-distance and the intrinsic and extrinsic 
3-curvature tensors of the 3-spaces E T ; 

b) the proper time of a time-like observer; 

c) the time-like and null 4-geodesics (they are relevant for the definition of 
the radial velocity of stars as shown in the IAU conventions of Rcf.[81] and 
in study of time delays [82], [80]); 

d) the red-shift and luminosity distance. In particular one finds that the 
recessional velocity of a star is proportional to its luminosity distance from 
the Earth at least for small distances. This is in accord with the Hubble old 
red-shift-distance relation which is formalized in the Hubble law (velocity- 
distance relation) when the standard cosmological model is used (see for 
instance Ref.[46] on these topics). These results have a dependence on the 
non-local York time, which could play a role in giving a different interpre- 
tation of the data from super-novae, which are used as a support for dark 
energy [1]. 



The properties of HPM transverse electro-magnetic fields have still to be explored. 
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Finally, in Subsection IIIB of the second paper in Refs.[40] it is shown 
that this HPM linearization can be interpreted as the first term of a HPM 
expansion in powers of the Newton constant G in the family of 3-orthogonal 
gauges. This expansion has still to be studied. In particular it will be useful 
to check whether in the HPM formulation there are phenomena (appearing at 
high orders in the standard PN expansions) like the hereditary tails starting 
from 1.5PN [0((^) 3 )] and the non-linear (Christodoulou) memory starting 
from 3PN (see Ref. [83] for a review) 36 . This would allow one to make a 
comparison with all the results of the PN expansions, in which today there is 
control on the GW solution and on the matter equations of motion till order 
3.5PN ) 7 )] (for binaries see the review in chapter 4 of Ref.[44]) and well 
established connections with numerical relativity (see the review in Ref. [84]) 
especially for the binary black hole problem (see the review in Ref. [85]). 



36 They imply that GW propagate not only on the flat light-cone but also inside it (i.e. 
with all possible speeds < v < c): there is an instantaneous wavefront followed by a 
tail traveling at lower speed (it arrives later and then fades away) and a persistent zero- 
frequency non-linear memory. 
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8 Post-Minkowskian Hamilton Equations for Particles, 
their Post-Newtonian Limit and Dark Matter as a 
Relativistic Inertial Effect 

The PM Hamilton equations and their PN limit in 3-orthogonal gauges for 
a system of N scalar particles of mass rrii and Grassmann-valued signs of 
energy % is discussed in this Section by using the results of the third paper in 
Ref.[40]. See Refs.[71, 78, 86, 87] for classical texts on the motion of particles 
in gravitational fields and Refs.[74, 88, 89] for more recent developments 37 . 

The treatment in the 3-orthogonal gauges of the PM Hamilton equations 
for the electro-magnetic field in the radiation gauge is given in the second 
paper of Ref.[40], while the PM Hamilton equations for perfect fluids are 
given in Ref.[41]. 

With only particles the PM approximation with the ultraviolet cutoff 
M implies Ki{r) = + McO(C), M {1) (r,a) - £. S 3 (a, Vi (r)) Vi 

7 JM= + 0(eiM (1)r V (T^) =Xi <5 3 (^r, 4 (r))^^|g+0(C 2 ).More- 

over one has %(r) = 0(C)- The notation d(r) = da ^ is used. 

One can make a equal time development of the retarded kernel in Eq.(54) 
like in Ref.[27] for the extraction of the Darwin potential from the Lienard- 
Wiechcrt solution (see Eqs. (5.1)-(5.21) of Ref.[27] with £ s p ± (<*) Vi (r) ^ 
J2uv 'Pbbuv(cr) ^1^2^ )- ^ n tms wa y one S ets tnc following expression of 
the HPM GW from point masses 



C b i uv V 1 -^^) 

J^U 1 / \2m 

m=l v ' 

+ o(c 2 ), 

d r d s \ . 1 /. d r d s \ d u d v 



+ 



A J A 



37 In this approach point particles arc considered as independent matter degrees of freedom 
with a Grassmann regularization of the self-energies to get well defined world-lines (see 
also Rcf. [89]): they are not considered as point-like singularities of solutions of Einstein's 
equations (the point of view of Ref. [86]). Solutions of this type have to be described with 
distributions and, as shown in Rcf. [90], the most general class of such solutions under 
mathematical control includes singularities simulating matter shells, but not either strings 
or particles. See also Ref. [91]. 
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d u d v 

(5 rv d s + 5 SV d r ) + (S ru d s + S su d r ) 



(59) 



with the retardation effects pushed to order 0(( 2 ). 

If the lapse and shift functions are rewritten in the form ri(i) = — 
d T 3 ^C(i), ^(i)(r) — ^(i)(r) + d r 3 ^(i)j to display their dependence on the iner- 
tial gauge variable non-local York time, it can be shown that the PM Hamil- 
ton equations for the particles imply the following form of the PM Grassmann 
regularized second-order equations of motion showing explicitly the equality 
of the inertial and gravitational masses of the particles 



m. 



def 



mJ 7 I(r\r] l (T)\r] k ^ l (T)) 



( 



+ 



(E* u w[ 



<9»(1)( M ) ^^(l)(r 



]-EEw 



a 77" 



+ 



9^ 



(T,r, i (r)) + (^(^(r)) 2 



9(jj 1) +2 (1) ) 



"* w E h w I 2 — — + V — Wt — ) 
-Ew ( 2 



9(r r w + 2 (1) ) 



2 — + 

TO 
l-^(r) 



(,.);.'(»)-' <;( r c (1) +2 (1) ) 
<9< 



02 1 a r , o \- dd r | q . 3 /C (1) 



drfi 
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drjfdrj: 



(60) 



The effective action-at-a-distance force Fi(r) contains 

a) the contribution of the lapse function , which generalizes the Newton 
force; 

b) the contribution of the shift functions «(i)( r ), which gives the gravito- 
magnetic effects; 

c) the retarded contribution of HPM GW, described by the functions 
of Eq.(59); 

d) the contribution of the volume element <f>^ ((f) — 1 + 6 0(i) + 0(C 2 )), 
always summed to the HPM GW, giving forces of Newton type; 

e) the contribution of the inertial gauge variable (the non-local York time) 

3/C (l) = A" 3X (1)- 

In the electro-magnetic case in SR [28] the regularized coupled second- 
order equations of motion of the particles obtained by using the Lienard- 
Wiechert solutions for the electro-magnetic field are independent by the type 
of Green function (retarded or advanced or symmetric) used. The electro- 
magnetic retardation effects, killed by the Grassmann regularization, are con- 
nected with QED radiative corrections to the one-photon exchange diagram. 
This is not strictly true in the gravitational case. The effect of retardation 
is not killed by the Grasmann regularization but only pushed to 0(( 2 ): at 
this order it should give extra contributions to the second-order equations of 
motion. This shows that our semi-classical approximation, obtained with our 
Grassmann regularization, of a unspecified "quantum gravity" theory does 
not take into account only a " one-graviton exchange diagram" : in the spin 
2 case there is an extra retardation effect showing up only at higher HPM 
orders 38 . 



8.1 The Center-of-Mass Problem in General Relativity 
and in the HPM Linearization. 

As said in Section 5, the 3-universe is described in a non inertial rest frame 
with non-Euclidean 3-spaces S T tending to Euclidean inertial ones ZV(oc) 



38 In the electro- magnetic case the Grassmann regularization implies QiT)\{j — |<x|) = 
QiVi( T ) an d equations of motion of the type rjV(r) = Q; ... with Q| = 0. In the gravita- 
tional case the equations of motion are of the type r\i r\\(r) = r]i ... with r)f = 0, but the 
Grassmann regularized retardation in Eq.(54) gives Eq.(59) only at the lowest order in f 
and has contributions of every order 0(( k ). 
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at spatial infinity. Both matter and gravitational degrees of freedom live 
inside S T and their internal 3-center of mass is eliminated by the rest-frame 
condition P ADM ~ (implied by the absence of super-translations) if also 
the condition K ADM = J T ADM rts is added like in SR. The 3-universe may 
be described as an external decoupled center of mass carrying a pole-dipole 
structure: Eadm is the invariant mass and J ADM the rest spin. As in SR 
the condition K r ADM k, selects the Fokker=Prycc center of inertia as the 
natural time-like observer origin of the radar coordinates: it follows a non- 
geodetic straight world-line like the asymptotic inertial observers existing in 
these space-times. 

This is a way out from the the problem of the center of mass in general 
relativity and of its world-line, a still open problem in generic space-times as 
can be seen from Refs. [58, 92] (and Ref.[74] for the PN approach). Usually, 
by means of some supplementary condition, the center of mass is associated 
to the monopole of a multipolar expansion of the energy-momentum of a 
small body (see Ref.[26] for the special relativistic case). 

In SR the elimination of the internal 3-center of mass leads to describe the 
dynamics inside S T only in terms of relative variables (see Eqs.(15) in the 
case of particles) . However relative variables do not exist in the non-Euclidean 
3-spaccs of curved space-times, where flat objects like rij(r) = rj^r) — Vj{ T ) 
have to be replaced with a quantity proportional to the tangent vector to 
the space-like 3-geodesics joining the two particles in the non-Euclidean 3- 
space E T (see Ref. [93] for an implementation of this idea). Quantities like 
rfj(r) have to be replaced with the Synge world function [82, 87, 89] 39 . This 
problem is another reason why extended objects tend to be replaced with 
point-like multipoles, which, however, do not span a canonical basis of phase 
space (see Refs. [26] for SR). 

However, at the level of the HPM approximation one can introduce relative 
variables for the particles, like the SR ones of Eq.(15), defined as 3- vectors in 
the asymptotic inertial rest frame ZV(oc) by putting tj^t) — V(o)i( T )+ r l{i)i(' T ) 
and Ki(r) = k {o)i (t) + k {1)i (t) with tj (o)4 (t), k {o)i (t) = 0(C). This allows 
one to define HPM collective and relative canonical variables for the particles, 
with the collective variables eliminated by the conditions PJ flM « and 
K r ADM « (at the lowest order they become the SR conditions). 

In the case of two particles (with total and reduced masses M — mi + m 2 
and n = ™±p.) one puts ^(r) = t] 12 (t) + ^ p 12 (r), t) 2 {t) = ij 12 (t) - 
ffeW, KiW = ^K 12 (T)+ir 12 (T), k 2 (t) = ^Ki2(t)-7Ti2(t) and goes 
to the new canonical basis ^^(t) ~ mi T1 ^ T ^ m2 ^ ^ p 12 (r) — »7 1 (r) — J7 2 ( T )i 
«i 2 (r) = Ki(r) + « 2 (r), tt 12 (t) = 



39 It is a bi-tensor, i.e. a scalar in both the points ■q i (r) and r)^{r), defined in terms of the 
space-like geodesic connecting them in E T . See Eq.(3.13) of the third paper in Ref. [40]. 



68 Luca Lusanna 

It can be shown that the conditions P T ADM ~ and K r ADM w imply 



■>1i{t) ~ - A (o) (t)) p (o)12 (r) + ^P ( i)i 2 (t) + f {1) (T)[rel.var.,GW], 
V2( T ) ~ -{jf + A (o)( T )) P(o)i 2 ( T ) - ^-P(i)i2( T ) + f (1) (T)[rel.var.,GW], 

m M c2 + - ^ M ^ + ^i)i 2 w 

A (o)( T ) = / = , = , 

\j m \ C 2 + «fl)l 2 ( T ) + \/ m 2 C 2 + W (1)1 2 ( T ) 

(61) 

for some function [rel.var., GW](t) rts V(i)i2( T ) depending on the relative 
variables and the HPM GW of Eq.(59) in absence of incoming radiation. Then 
the equations of motion (60) imply 



^P(o)i 2 ( T ) = f i(t\t)( )i(t)\V(o)2{t)) - -j^F^(T\r 1{o)2 (T)\r 1{on (T)), 

(62) 

for the relative configurational variable. The collective configurational vari- 
able has T7( )i 2 (t) ~ — j4( q )(t) P( )i 2 ( t ) a t the lowest order, while at the first 
order there is an equation of motion equivalent to ^[ 1 ) 12 ( T ) ~ £pz /(i) [rel.var., GW] 
(r). 



8.2 The Post- Newtonian Expansion at all Orders in 
the Slow Motion Limit. 



If all the particles are contained in a compact set of radius l c , one can add a 
slow motion condition in the form ^Je = ^ « \J -y 24 -, i = 1, --iV (i? mi = 2 ^j"* 
is the gravitational radius of particle i) with l c > R M and A >> l c (see the 
Introduction). In this case one can do the PN expansion of Eqs.(60). 

After having put r = ct, one makes the following change of notation 



Vi(r) = Vi(t), Vi(t) = 

Vi(t) 



dt 



a-i (<) 



dt 2 ' 



(63) 
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For the non-local York time one uses the notation 3 IC(i) (t, <r) 
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= 3 /C (1) (r,<T). 



Then one studies the PN expansion of the equations of motion (60) with 
the result (kPN means of order 0(c~ 2k )) 



d 2 ffj{t) r d x ^ mj 1 ^[(*) /o2. 3^ , 

+ 2 E *w + E *w ) (*. *(*)) 

+ ^(ipn) (*) + (Wp/ier PN orders). (64) 
At the lowest order one finds the standard Newton gravitational force 

Fi(Newton)(t) = -TUi G g-^r {t)-r) j (t) \ ' 

The unexpected result is a 0.5PN force term containing all the dependence 
upon the non-local York time. The (arbitrary in these gauges) double rate 
of change in time of the trace of the extrinsic curvature creates a 0.5 PN 
damping (or anti-damping since the sign of the inertial gauge variable 3 IC(i) 
of Eq.(57) is not fixed) effect on the motion of particles. This is a inertial 
effect (hidden in the lapse function) not existing in Newton theory where the 
Euclidean 3-space is absolute. 

Then there are all the other kPN terms with k = 1, 1.5, 2,... Since these 
results have been obtained without introducing ad hoc Lagrangians for the 
particles, are not in the harmonic gauge and do not contain terms of order 
0(£ 2 ) and higher, it is not possible to make a comparison with the standard 
PN expansion (whose terms are known till the order 3.5PN [44]). Therefore 
only the 1PN and 0.5PN terms will be considered in the next two Subsections. 



8.3 The HPM Binaries at the 1PN Order 

Since in the next Subsection the 0.5PN term depending on the non-local 
York time will be connected with dark matter at the level of galaxies and 
clusters of galaxies and since there is no convincing evidence of dark matter 
in the Solar System and near the galactic plane of the Milky Way [96], it is 
reasonable to assume 3 ^C(i) (t, a) = (r, er) ss near a star with planets 

and near a binary. 

In the description of Subsection 8.1 of the 1PN two-body problem, which 
is relevant for the treatment of binary systems 40 as shown in Chapter VI 

40 For binaries one assumes ^ fa -5^ << 1, where l c \r\ with f(t) being the relative 
separation after the decoupling of the center of mass. Often one considers the case mi Ri 
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of Refs.[44] based on Ref.[47, 94, 95], it can be shown that the relative mo- 
mentum in the rest frame has the 1PN expression 7t 12 (t) = 7T( 1 ) 12 (t) w 



, m 1 +m' 2 ^( re ;)(o)i2(*) ^2 
1 + 2 MS I E J 



where v 



(rei)(o)12 



(*) 



d P(o)12(*) 



IS 



the velocity of the lowest order P( )i 2 (r) of the relative variable. 

If one ignores the York time and considers only positive energy particles 
(771,772 !-> +1), the 1PN equations of motion for the relative variable of the 
binary implied by Eqs.(62) and 1PN expression of the weak ADM energy 
Eadm and of the rest spin J A S DM (being determined by the boundary condi- 
tions they are constants of the motion implying planar motion in the plane 
orthogonal to the rest spin ) can be shown to be 



dV(rel)(o)12(t) 

~dt 



-GM 



P(o)12(*) 



1+,1+3 M- 



M \ u (rel)(o) 



\p(o)12(tW 
3/J ( V(rel)(o)12{t) ■ P( o )i 2 (0) \ 
^ |P(o)12(i)l I 



o)l 2 (*) 



2M 



+ 



GM ( 2/j\ u v {rel)(o)12 (t) ■ p {o)12 (t)) 

~ I^WFV 4 - mJ-m)( 0)12 ( t) • 

E ADM{1PN) = £ mi c 2 + , (i < e0(o)12 (*) [l + f 



GM 



|p(o)l 2 WI 



1 + 



V frel) (0)12(0 



+ 



Trs 

J ADM(IPN) 



V ( V(rel)(o)12(t) P(o)\l{t) . 2 \1\ 

m { c ' i P(0)12 (i)r ^ 

(P( )12(*)^ eO(o)12W P(o)12W V (re0(o)12 
mf + m| ( V(rel)(o)12(t) ]2 

2M 3 1 c j 



(65) 



Our 1PN equations (65) in the 3-orthogonal gauges coincide with Eqs. 
(2.5), (2.13) and (2.14) of the first paper in Ref.[47] (without G 2 terms since 
they are 0*(£ 2 )), which are obtained in the family of harmonic gauges starting 
from an ad hoc 1PN Lagrangian for the relative motion of two test particles 
(first derived by Infeld and Plebanski [86] ) 41 . These equations are the start- 
ing point for studying the post-Keplerian parameters of the binaries, which, 

m2. See chapter 4 of Rcf. [44] for a review of the emission of GW's from circular and elliptic 
Keplerian orbits and of the induced inspiral phase. 

41 This is also the starting point of the effective one body description of the two-body 
problem of Rcfs. [88]. 
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together with the Roemer, Einstein and Shapiro time delays (both near Earth 
and near the binary) in light propagation, allow one to fit the experimental 
data from the binaries (see the second paper in Ref.[47] and Chapter VI 
of Ref. [44] ) . Therefore these results are reproduced also in our 3-orthogonal 
gauge with 3 /C(i)(t, <t) = 0. 



8-4 From the Three Signatures for Dark Matter 
Reinterpreted as Relativistic Inertial Effects 
Induced by the York Time to the Need of a PM 
ICRS 

To study the effects induced by the 0.5PN velocity-dependent (friction or anti- 
friction) force term in Eqs.(64), depending on the inertial gauge variable non- 
local York time 3 IC( 1 )(t, a) — 3 K^(t, er) w ^ F^(t,<t) with F^(t,<t) 
arbitrary numerical function, it is convenient to rewrite such equations in the 
form 



d_ 

dt 



+ c dt ^ ) -^T ) = ~ G dm E * \*M-Lm + 



d 



+ 0(C 2 ), (66) 



because the damping or anti-damping factors in Eq.(64) are 7t(i,Tjj(t)) = 
^% 1) (t,f} i (t))a.ndfj i (t) = 0(C). 

As a consequence the velocity-dependent force can be reinterpreted as the 
introduction of an effective (time-, velocity- and position-dependent) inertial 
mass term for the kinetic energy of each particle: 



m t ^mi(l + ^ ^ 3 ^)(i)(t, r/i(t)) ^) = m t + (A m)i(t, fj^t)), 

(67) 

in each instantaneous 3-space. Instead in the Newton potential there are 
the gravitational masses of the particles, equal to the inertial ones in the 4- 
dimensional space-time due to the equivalence principle. Therefore the effect 
is due to a modification of the effective inertial mass in each non- Euclidean 3- 
space depending on its shape as a 3-sub-manifold of space-time: it is the equal- 
ity of the inertial and gravitational masses of Newtonian gravity to be violated}. 
In Galilei space-time the Euclidean 3-space is an absolute time-independent 
notion like Newtonian time: the non-relativistic non-inertial frames live in 
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this absolute 3-space differently from what happens in SR and GR, where 
they are (in general non-Euclidean) 3-sub-manifolds of the space-time. 

Eqs. (64), (66) and (67) can be applied to the three main signatures of the 
existence of dark matter in the observed masses of galaxies and clusters of 
galaxies, where the 1PN forces are not important, namely the virial theorem 
[49] , the weak gravitational lensing [50] , [49] and the rotation curves of spiral 
galaxies (see Ref. [48] for a review) , to give a reinterpretation of dark matter 
as a relativistic inertial effect. 

A) Masses of clusters of galaxies from the virial theorem. For a bound sys- 
tem of N particles of mass m (N equal mass galaxies) at equilibrium, the virial 
theorem relates the average kinetic energy < Ekin > in the system to the av- 
erage potential energy < U po t > in the system: < Ekin >= — \ < U po t > 
assuming Newton gravity. For the average kinetic energy of a galaxy in the 
cluster one takes < Ekin >~ \m < v 2 >, where < v 2 > is the average of the 
square of the radial velocity of single galaxies with respect to the center of 
the cluster (measured with Doppler shift methods; the velocity distribution 
is assumed isotropic) . The average potential energy of the galaxy is assumed 
of the form < U po t >~ — G 2 ^-, where M = Nm is the total mass of the 
cluster and 1Z — a R is a " effective radius" depending on the cluster size R 
(the angular diameter of the cluster and its distance from Earth are needed 
to find R) and on the mass distribution on the cluster (usually a w 1/2). 
Then the virial theorem implies M rts || < v 2 >. It turns out that this mass 
M of the cluster is usually at least an order of magnitude bigger that the 
baryonic matter of the cluster M\, ar — Nm (spectroscopically determined) . 
By applying Eqs. (64) to the equilibrium condition for a self-gravitating sys- 
tem, i.e. 2pJ2i m i I Vi(t) | 2 = with rm = m, one gets £\ m i w iW ~ 
G Hi>j \ns)-f) j {t)\ - 1 c Ei m i (*7i(*)'«i(*)) 7i(*,»ji(*)) = with m t = m 3 = 
m. Therefore one can write (U pot ) = -± J2 t>J ^fi™^ « G (with 
K = R/2) and \ m (v 2 ) = -\ (U pot ) + f c ((fj ■ v^j j(t, fj)) with the notation 

((fj-v^J l(t,fj)) = ^ J2i (jii(t) ■ liitiViit)) (it contains the non-local 

York time). Therefore for the measured mass M (the effective inertial mass 
in 3-space) one has 

M = ^ < 1,2 >= M bar + ^((^- V ) V)) d =M bar + M DM , 

(68) 

and one sees that the non-local York time can give rise to a dark matter 
contribution Mdm = M — Mb ar - 

B) Masses of galaxies or clusters of galaxies from weak gravitational lens- 
ing. Usually one considers a galaxy (or a cluster of galaxies) of big mass M 
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behind which a distant, bright object (often a galaxy) is located. The light 
from the distant object is bent by the massive one (the lens) and arrives 
on the Earth deflected from the original propagation direction. As shown 
in Ref.[50] one has to evaluate Einstein deflection of light, emitted by a 
source S at distance ds from the observer O on the Earth, generated by 
the big mass at a distance do from the observer O. The mass M, at dis- 
tance dDS from the source S, is considered as a point-like mass generating 
a 4- metric of the Schwarzschild type (Schwarzschild lens). The ray of light 
is assumed to propagate in Minkowski space-time till near M, to be de- 
flected by an angle a by the local gravitational field of M and then to prop- 
agate in Minkowski space-time till the observer O. The distances ds, do, 
dos-, & re evaluated by the observer O at some reference time in some nearly- 
inertial Minkowski frame with nearly Euclidean 3-spaces (in the Euclidean 
case dns = ds — do). If £ = 9 do is the impact parameter of the ray of light 
at M and if £ >> R s — 2< ^ M (the gravitational radius), Einstein's deflection 
angle is a — = i f 2 1 ^ and the so-called Einstein radius (or characteristic 

angle) is a Q = ^2 R s = A measurement of the deflection 

angle and of the three distances allows to get a value for the mass M of the 
lens, which usually turns out to be much larger of its mass inferred from the 
luminosity of the lens. For the calculation of the deflection angle one consid- 
ers the propagation of ray of light in a stationary 4-metric of the BCRS type 
and uses a version of the Fermat principle containing an effective index of re- 
fraction n. One has n — A g TT = e [1 — ^nr — 2 d T 3 /C] in the PM approximation. 

Since one has ^ = , the definition 2d T 3 IC (1) d = - g J?$ L leads to 

an Einstein deflection angle 

a= ±GM wM M t/M 6ar + M DM . (69) 
c s 

Therefore also in this case the measured mass M is the sum of a baryonic 
mass Mb ar and of a dark matter mass Mjjm induced by the non-local York 
time at the location of the lens. 

C) Masses of spiral galaxy masses from their rotation curves. In this case 
one considers a two-body problem (a point-like galaxy and a body circulating 
around it) described in terms of an internal center of mass ff^it) ~ V(i)i2(^) 
(rj( ) 12 (t) = is the origin of the 3-coordinates) and a relative variable Pi 2 (t)- 
Then the sum and difference of Eqs.(64) imply the equations of motion for 
V(i)i2(t) an d Pi2{t)- While the first equation implies a small motion of the 
overall system, the second one has the form 

d2 P(i)i 2 W o ~ P \ 2 {t) 1 d~p\ 2 {t) 

~ GM i- u\n 11 — 7+(*>Pi2W)«W;> 



dt 2 |pi2(i)| 3 c dt 



74 



Luca Lusanna 



7+ (i, pi 2 (t), V(t)) = -jj 71 (*, ^ Pl2(t), »(*)) + 72(*, --^ Pl2(t), »(*)), 

(70) 

where 7, are the damping or anti-damping factors defined after Eq.(66). 
Eq.(70) gives the two-body Kepler problem with an extra perturbative 
force. Without it a Keplerian solution with circular trajectory such that 
I Pi2(t) \= R = const, implies that the Keplerian velocity v a (t) = v n(t) 

has the modulus vanishing at large distances, v = \J~^jf- -^r^oo 0. Instead 
the rotation curves of spiral galaxies imply that the relative 3-velocity goes 
to constant for large R, i.e. v = G ( M t><^+ AM ( r )) — const. (Mb ar is 
the spectroscopically determined baryon mass), so that the extra required 
term AM{r) is interpreted as the mass Mdm of a dark matter halo. 

The presence of the extra force term implies that the velocity must be 
written as v(t) = v (t) + vi(t) with vi(t) a first order perturbative correction 
satisfying —j^— = — n(t) 7+(i, p\2{t),v r {t)). Therefore at the first order in 
the perturbation one gets v 2 (t) = v 2 ^ 1 — | h(t) ■ dt\ h(t\) 7+(ti, Pi2(ti), 

v D (ti))^j. Therefore, after having taken a mean value over a period T (the 

time dependence of the mass of a galaxy is not known) the effective mass of 
the two-body system is 



M eff = ^-^=M^l-(^h(t)- J dt 1 h{t 1 ) 1+ (t 1 ,p 12 {t 1 ),v {t 1 ))j^) = 

= M bar + M DM . (71) 

with a AM(r) — Mdm function only of the mean value of the total time 
derivative of the non-local 3 /C(i) to be fitted to the experimental data. 

Therefore, the existence of the inertial gauge variable York time, a prop- 
erty of the non-Euclidean 3-spaces as 3-sub-manifolds of Einstein space-times 
(connected only to the general relativistic remnant of the gauge freedom in 
clock synchronization, independently from cosmological assumptions) implies 
the possibility of describing part (or maybe all) dark matter as a relativistic 
inertial effect in Einstein gravity without alternative explanations using: 

1) the non-relativistic MOND approach [97] (where one modifies Newton 
equations); 

2) modified gravity theories like the f(R) ones (see for instance Refs.[98]; 
here one gets a modification of the Newton potential); 



3) the assumption of the existence of WIMP particles [99]. 
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Let us also remark that the 0.5PN effect has origin in the lapse function 
and not in the shift one, as in the gravito-magnetic elimination of dark matter 
proposed in Ref.[100]. 

The open problem with this explanation of dark matter is the determina- 
tion of the non-local York time from the data on dark matter. From what 
is known about dark matter in the Solar System and inside the Milky Way 
near the galactic plane, it seems that 3 /C(i) (t, er) is negligible near the stars 
inside a galaxy. Instead the non-local York time (or better a mean value in 
time of its total time derivative) should be relevant around the galaxies and 
the clusters of galaxies, where there are big concentrations of mass and well 
defined signatures of dark matter. Instead there is no indication on its value 
in the voids existing among the clusters of galaxies. 

Therefore the known data on dark matter do not allow one to get an 
experimental determination of the York time 3 K^(t, er) = A 3 /C(i)(t, er), 
because to do it one needs to know the non-local York time on all the 3- 
universe at a given r. 

Since, as said in the Introduction, at the experimental level the descrip- 
tion of matter is intrinsically coordinate- dependent, namely is connected with 
the conventions used by physicists, engineers and astronomers for the mod- 
eling of space-time, one has to choose a gauge (i.e. a 4-coordinate system) in 
non-modified Einstein gravity which is in agreement with the observational 
conventions in astronomy. This way out from the gauge problem in GR re- 
quires a choice of 3-coordinates on the instantaneous 3-spaces identified by 
a choice of time and by a clock synchronization convention, i.e. a fixation 
of the York time 3 K( 1 - ) (t,ct). The convention resulting by one set of such 
choices would give a PM extension of ICRS, with BCRS being its quasi- 
Minkowskian approximation for the Solar System. Since the existing ICRS 
[3, 5] has diagonal 3-metric, 3-orthogonal gauges are a convenient choice. 

The real problem is the extraction of an indication of which kind of 
function of time and 3-coordinates to use for the York time 3 Kn\(r, a) 
from astrophysical data different from the ones giving information about 
dark matter. Once one would have a phenomenological parametrization of 
the York time, then the data on dark matter would put restrictions on 
the induced phenomenological parametrization of the non-local York time 
3 /C(!)(t, er) = 3 if( 1 )(r, er). As it will be delineated in the final Section, to 
implement this program one has to look at the astrophysical data on dark en- 
ergy after having succeeded to interpret also it as a relativistic inertial effect 
in suitable cosmological space-times in which one can induce the distinction 
between inertial and tidal degrees of freedom of the gravitational field from 
the previously discussed Hamiltonian framework. 
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9 Dark Energy and Other Open Problems 

This Lecture contains a full review of an approach to SR and to asymptoti- 
cally Minkowskian classical canonical Einstein GR based on a description of 
global non-inertial frames centered on a time- like observer which is suggested 
by relativistic metrology. The gauge freedom in clock synchronization, which 
does not exist in Galilei space-time (Newton time and Euclidean 3-spaces are 
absolute) and is not restricted in Minkowski space-time (it spans the class 
of the admissible 3+1 splittings of this absolute space-time), is restricted in 
GR to the gauge freedom connected with the inertial gauge variable 3 K , the 
York time, which determines the shape of the instantaneous non-Euclidean 
3-spaces as 3-sub-manifolds of the space-time. 

The study of canonical ADM tetrad gravity in asymptotically Minkowskian 
space-times without super-translations (so that they admit an asymptotic 
ADM Poincare algebra at spatial infinity) in the York canonical basis al- 
lowed one to disentangle the tidal degrees of freedom of the gravitational 
field from the inertial gauge ones (they include the York time), to find the 
family of non-harmonic 3-orthogonal Schwinger time gauges and to define 
a HPM linearization in them. The main properties of these non-harmonic 
gauges are that only the HPM-GW (but not the lapse and shift functions) 
are retarded quantities with a no-incoming radiation condition and that one 
can naturally find which quantities depend upon the York time. 

Relativistic particle mechanics, coupled to the electro-magnetic field in 
the radiation gauge, has been studied both in SR and GR with a suitable 
Grassmann regularization of the self-energies so to get well defined equations 
of motion. 

In SR, after a clarification of the problem of the relativistic center of mass 
and the definition of inertial and non-inertial rest frames of isolated systems, 
it was possible to develop a formulation, the parametrized Minkowski the- 
ories, in which the transitions among global non-inertial frames are gauge 
transformations. Then isolated systems were described in the rest-frame in- 
stant form of dynamics and the structure of their Poincare generators and 
of their relative variables in the instantaneous Wigner 3-spaces was clarified. 
With this approach it was possible to give a new formulation of the micro- 
canonical ensemble in relativistic kinetic theory and to develop a formulation 
of relativistic quantum mechanics and relativistic entanglement taking into 
account the known results about relativistic bound states and the spatial non- 
separability and non-locality induced by the Lorentz signature of Minkowski 
space-time. 

In GR it was possible to derive regularized equations of motion of the 
particles in the non-inertial rest frame and to study their PM limit in the 
HPM linearization in the 3-orthogonal gauges and the emission of HPM GW 
(with the energy balance under control even in absence of self- forces) . Then 
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the PN limit of these PM equations allows one to recover the known 1PN 
results of harmonic gauges. The more surprising result is that in the PN 
expansion of the PM equations of motion there is a 0.5PN term in the forces 
depending upon the York time. This opens the possibility to describe dark 
matter as a relativistic inertial effect implying that the effective inertial mass 
of particles in the 3-spaces is bigger of the gravitational mass because it 
depends on the York time (i.e. on the shape of the 3-space as a 3-sub-manifold 
of the space-time: this is impossible in Newton gravity in Galilei space-time 
and leads to a violation of the Newtonian equivalence principle). 

The proposed solution to the gauge problem in GR based on the conven- 
tions of relativistic metrology for ICRS and the results of the last Section on 
the re-interpretation of dark matter as a relativistic inertial effect arising as 
a consequence of a convention on the York time in an extended PM ICRS 
push toward the necessity of similar re-interpretation also of dark energy in 
cosmology [1, 101, 102, 103]. As it has been shown, the identification of the 
tidal and inertial degrees of freedom of the gravitational field can be reformu- 
lated in the framework of the non-Hamiltonian first-order ADM equations by 
means of the replacement of the Hamiltonian momenta with the expansion 
and the shear of the Eulerian observers associated with the 3+1 splitting of 
the space-time. Therefore this identification can also be applied to the cos- 
mological space-times which do not admit a Hamiltonian formulation: also 
in them the identification of the instantaneous 3-spaces E T , now labeled by 
a cosmic time, requires a conventional choice of clock synchronization, i.e. a 
convention on the York time 3 K defining the shape of the 3-spaces as 3-sub- 
manifolds of the space-time, and of 3-coordinates (the 3-orthogonal ones are 
acceptable also in cosmology). 

In the standard ylCDM cosmological model the class of cosmological so- 
lutions of Einstein equations is restricted to Friedmann-Robertson- Walker 
(FRW) space-times with nearly Euclidean 3-spaces (i.e. with a small internal 
3-curvature) . In them the Killing symmetries connected with homogeneity 
and isotropy imply (r is the cosmic time, a(r) the scale factor) that the York 
time is no more a gauge variable but coincides with the Hubble constant: 
3 K (t) = — = —H(t). However at the first order in cosmological pertur- 
bations (see Rcf. [104] for a review) one has 3 K = —H + 3 K/u with 3 K^ 
being again an inertial gauge variable to be fixed with a metrological conven- 
tion. Therefore the York time has a central role also in cosmology and one 
needs to know the dependence on it of the main quantities, like the red-shift 
and the luminosity distance from supernovae, which require the introduction 
of the notion of dark energy to explain the 3-universe and its accelerated 
expansion in the framework of the standard ylCDM cosmological model. 

Instead in inhomogeneous space-times without Killing symmetries like the 
Szekeres ones [105, 106] the York time remains an arbitrary inertial gauge 
variable. Therefore the main open problem of the present approach is to 
see whether it is possible to find a 3-orthogonal gauge in a inhomogeneous 
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Einstein space-time (at least in a PM approximation) in which the convention 
on the inertial gauge variable York time allows one to accomplish the following 
two tasks simultaneously: a) to eliminate both dark matter and dark energy 
through the choice of a 4-coordinate system (suggested by astrophysical data) 
to be used in a consistent PM reformulation of ICRS and b) to save the main 
good properties of the standard ACDM cosmological model due to the inertial 
and dynamical properties of the space-time. As matter one will take the dust, 
whose description in the York canonical basis is given in Ref.[41]. 

Also in the back-reaction approach [107, 108] to cosmology, according to 
which dark energy is a byproduct of the non-linearities of GR when one con- 
siders spatial averages of 3-scalar quantities in the 3-spaces on large scales 
to get a cosmological description of the universe taking into account its ob- 
served inhomogeneity, one gets that the spatial average of the product of 
the lapse function and of the York time (a 3-scalar gauge variable) gives the 
effective Hubble constant. Since this approach starts from the Hamiltonian 
description of an asymptotically flat space-time and since all the canonical 
variables in the York canonical basis, except the angles 6 l , are 3-scalars, the 
formalism presented in this Lecture will allow to study the spatial average 
of nearly all the Hamilton equations and not only of the super-Hamiltonian 
constraint and of the Hamilton equation for the York time as in the existing 
formulation of the approach. This will be done by using the perfect fluids of 
Ref.[41] as matter. 

Also the recent point of view of Rcf.[109], taking into account the relevance 
of the voids among the clusters of galaxies, has to be reformulated in terms 
of the York time. 

Finally one should find the dependence upon he York time of the Landau- 
Lifschitz energy-momentum pseudo-tensor and re-express it as the effective 
energy-momentum tensor of a viscous pseudo-fluid. One will have to check 
whether for certain choices of the York time the resulting effective equation of 
state of the fluid has negative pressure, realizing also in this way a simulation 
of dark energy. 

Other open problems in GR under investigation are: 

A) Find the second order of the HPM expansion to see whether in PM 
space-times there is the emergence of hereditary terms [44, 83] like the ones 
present in harmonic gauges. 

B) Study the PM equations of motion of the transverse electro-magnetic 
field trying to find Lienard-Wiechert-type solutions in GR. Study astrophys- 
ical problems where the electro-magnetic field is relevant. 

C) Find the expression in the York canonical basis of the Weyl scalars of 
the Newman-Penrose approach [68] and then of the four Weyl eigenvalues, 
which are tetrad-independent 4-scalar invariants of the gravitational field. 
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Is it possible to find a canonical transformation replacing the 3-scalar tidal 
variables with four 4-scalar functions of the Weyl eigenvalues? Are Weyl 
eigenvalues Dirac observables? 

D) Try to make a multi-temporal quantization (see Refs,[34, 65]) of the 
linearized HPM theory over the asymptotic Minkowski space-time, in which, 
after a Shanmugadhasan canonical transformation to a new York canonical 
basis adapted to all the constraints, only the tidal variables are quantized 
but not the inertial gauge ones. After this type of quantization, in which 
the lapse and shift functions remain c-numbers, the space-time would still 
be a classical 4-manifold: only the two eigenvalues of the 3-metric describing 
GW are quantized and therefore only 3-metric properties like 3-distances, 
3-areas, 3-volumes become quantum properties. After having re-expressed 
the Ashtekar variables [110] for asymptotically Minkowskian space-times (see 
Appendix B of Ref. [4] ) in this final York canonical basis it will be possible to 
compare the outcomes of this new type of quantization with loop quantum 
gravity. 
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